(* Content-type: application/vnd.wolfram.mathematica *)

(*** Wolfram Notebook File ***)
(* http://www.wolfram.com/nb *)

(* CreatedBy='Mathematica 12.1' *)

(*CacheID: 234*)
(* Internal cache information:
NotebookFileLineBreakTest
NotebookFileLineBreakTest
NotebookDataPosition[       158,          7]
NotebookDataLength[    132941,       2410]
NotebookOptionsPosition[    126741,       2330]
NotebookOutlinePosition[    127084,       2345]
CellTagsIndexPosition[    127041,       2342]
WindowFrame->Normal*)

(* Beginning of Notebook Content *)
Notebook[{

Cell[CellGroupData[{
Cell[BoxData[
 RowBox[{
  RowBox[{"(*", 
   RowBox[{"Horan", " ", "and", " ", 
    RowBox[{"Shortle", ".", " ", "The"}], " ", "Song", " ", "Remains", " ", 
    "Not", " ", "the", " ", 
    RowBox[{"Same", ":", " ", 
     RowBox[{
     "Correlated", " ", "Intercept", " ", "and", " ", "Slope", " ", 
      "Uncertainties", " ", "Matter", " ", "to", " ", "Prices", " ", "vs", 
      " ", 
      RowBox[{"Quantities", ".", "  ", "JAERE", ".", " ", "2020."}], "  ", 
      "\[IndentingNewLine]", "Numerical", " ", "Example", " ", "1"}], ":", 
     " ", 
     RowBox[{"Quadratic", " ", "Restricted", " ", "Profits"}]}]}], "*)"}], 
  "\[IndentingNewLine]", 
  RowBox[{"Quit", "[", "]"}]}]], "Input",
 CellChangeTimes->{{3.805469800932412*^9, 3.8054698152641153`*^9}, {
  3.8125440425487757`*^9, 3.812544138969357*^9}, {3.81254417918489*^9, 
  3.8125441890242352`*^9}},
 CellLabel->
  "In[106]:=",ExpressionUUID->"53974844-04cd-41a6-a1d0-a926fb069df5"],

Cell[BoxData["\<\"Quadratic Restricted Profits\"\>"], "Output",
 CellChangeTimes->{
  3.8054698272580643`*^9, 3.8054701979282804`*^9, 3.805470236183175*^9, 
   3.805470269407214*^9, {3.805470309187398*^9, 3.805470337177867*^9}, 
   3.805470370768553*^9, 3.8054706565692344`*^9, 3.805471934963024*^9, 
   3.8054723800000467`*^9, 3.8054724719604907`*^9, 3.805472689436198*^9, 
   3.805473240510559*^9, {3.8054732804638376`*^9, 3.8054733095441556`*^9}, 
   3.805473859587235*^9, 3.8054739513185606`*^9, 3.805474018608631*^9, 
   3.8054745573527465`*^9, 3.805474837062269*^9, 3.80547543675031*^9, 
   3.805476465396245*^9, 3.805479179131485*^9, 3.805480107681384*^9, 
   3.805480421964707*^9, 3.8060612678780756`*^9, 3.8061435379992776`*^9, 
   3.8124798879765143`*^9},
 CellLabel->"Out[1]=",ExpressionUUID->"1e6e590f-bfc1-41dc-9797-22260b874c3d"]
}, Open  ]],

Cell[BoxData[
 RowBox[{
  RowBox[{"(*", 
   RowBox[{
    RowBox[{
    "Define", " ", "quadratic", " ", "restricted", " ", "profit", " ", 
     "function", " ", "and", " ", "damage", " ", 
     RowBox[{"function", ".", "  ", "Note"}], " ", "that", " ", "\[Pi]e", " ",
      "is", " ", "the", " ", "profit", " ", "function", " ", "restricted", 
     " ", "on", " ", "emisisons", " ", "rather", " ", "than", " ", "the", " ",
      "partial", " ", "derivative", " ", "of", " ", "profits", " ", "with", 
     " ", "respect", " ", "to", " ", 
     RowBox[{"emissions", ".", " ", "We"}], " ", 
     RowBox[{"can", "'"}], "t", " ", "define", " ", "\[Pi]", " ", "as", " ", 
     "profits", " ", "in", " ", "Mathematica", " ", "because", " ", 
     "Mathematica", " ", "only", " ", "recognizes", " ", "\[Pi]", " ", "as", 
     " ", "3.14"}], "..."}], "*)"}], "\[IndentingNewLine]", 
  RowBox[{
   RowBox[{
    RowBox[{"\[Pi]e", "=", 
     RowBox[{"\[Pi]e0", "+", 
      RowBox[{
       RowBox[{"(", 
        RowBox[{"\[Pi]e1", "+", "\[Eta]1"}], ")"}], "*", "e"}], "+", 
      RowBox[{"0.5", "*", 
       RowBox[{"(", 
        RowBox[{"\[Pi]e2", "*", "\[Eta]2"}], ")"}], "*", 
       RowBox[{"e", "^", "2"}]}]}]}], ";"}], "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"dam", "=", 
     RowBox[{"d1", "*", 
      RowBox[{"e", "^", "2"}]}]}], ";"}]}]}]], "Input",
 CellChangeTimes->{{3.805469843064266*^9, 3.805469852838595*^9}, {
  3.805469901966302*^9, 3.8054700044621105`*^9}, {3.8054725177738047`*^9, 
  3.8054725294881864`*^9}, {3.805473283339712*^9, 3.805473308126444*^9}, {
  3.8124748937980433`*^9, 3.812474920035285*^9}, {3.812476057561067*^9, 
  3.812476112310302*^9}},
 CellLabel->"In[1]:=",ExpressionUUID->"1fc03bfa-a798-4792-97da-c543b3a20ed5"],

Cell[BoxData[
 RowBox[{
  RowBox[{"(*", 
   RowBox[{
    RowBox[{
    "Compute", " ", "Expected", " ", "profits", " ", "assuming", " ", "a", 
     " ", "uniform", " ", "distribution", " ", "for", " ", "the", " ", "two", 
     " ", "independent"}], ",", " ", 
    RowBox[{
     RowBox[{"uncertain", " ", 
      RowBox[{"parameters", ":", " ", 
       RowBox[{
        RowBox[{"\[Eta]i", "~", "U"}], 
        RowBox[{"(", 
         RowBox[{"\[Eta]im", ",", "\[Eta]iM"}], ")"}], " ", "for", " ", 
        "i"}]}]}], "=", "1"}], ",", 
    RowBox[{
     RowBox[{
     "2.", " ", "The", " ", "distributions", " ", "will", " ", "be", " ", 
      "specified", " ", "such", " ", "that", " ", "E", 
      RowBox[{"{", "\[Eta]1", "}"}]}], " ", "=", " ", 
     RowBox[{
      RowBox[{"0", " ", "and", " ", "E", 
       RowBox[{"{", "\[Eta]2", "}"}]}], " ", "=", " ", "1"}]}]}], "*)"}], 
  "\[IndentingNewLine]", 
  RowBox[{"(*", 
   RowBox[{
   "Weighted", " ", "integral", " ", "of", " ", "profits", " ", "with", " ", 
    "respect", " ", "to", " ", "\[Eta]1"}], "*)"}], "\[IndentingNewLine]", 
  RowBox[{
   RowBox[{
    RowBox[{"E\[Pi]eint\[Eta]1", "=", 
     RowBox[{
      RowBox[{"Integrate", "[", 
       RowBox[{"\[Pi]e", ",", "\[Eta]1"}], "]"}], "/", 
      RowBox[{"(", 
       RowBox[{"\[Eta]1M", "-", "\[Eta]1m"}], ")"}]}]}], ";"}], 
   "\[IndentingNewLine]", "\[IndentingNewLine]", 
   RowBox[{"(*", 
    RowBox[{
     RowBox[{
     "Evaluate", " ", "integral", " ", "at", " ", "the", " ", "limits", " ", 
      "to", " ", "obtain", " ", "expectation", " ", "of", " ", "profits", " ",
       "with", " ", "respect", " ", "to", " ", "\[Eta]1"}], ",", " ", 
     RowBox[{"conditional", " ", "on", " ", "\[Eta]2"}]}], "*)"}], 
   "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"E\[Pi]e\[Eta]1", "=", 
     RowBox[{"Simplify", "[", 
      RowBox[{
       RowBox[{"(", 
        RowBox[{"E\[Pi]eint\[Eta]1", "/.", 
         RowBox[{"\[Eta]1", "\[Rule]", "\[Eta]1M"}]}], ")"}], "-", 
       RowBox[{"(", 
        RowBox[{"E\[Pi]eint\[Eta]1", "/.", 
         RowBox[{"\[Eta]1", "\[Rule]", "\[Eta]1m"}]}], ")"}]}], "]"}]}], 
    ";"}], "\[IndentingNewLine]", "\[IndentingNewLine]", 
   RowBox[{"(*", 
    RowBox[{
    "Weighted", " ", "integral", " ", "of", " ", "E\[Pi]e\[Eta]1", " ", 
     "with", " ", "respect", " ", "to", " ", "\[Eta]2"}], "*)"}], 
   "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"E\[Pi]eint\[Eta]2", "=", 
     RowBox[{
      RowBox[{"Integrate", "[", 
       RowBox[{"E\[Pi]e\[Eta]1", ",", "\[Eta]2"}], "]"}], "/", 
      RowBox[{"(", 
       RowBox[{"\[Eta]2M", "-", "\[Eta]2m"}], ")"}]}]}], ";"}], 
   "\[IndentingNewLine]", "\[IndentingNewLine]", 
   RowBox[{"(*", 
    RowBox[{"Overall", " ", "expected", " ", "profits"}], "*)"}], 
   "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"E\[Pi]e", "=", 
     RowBox[{"Simplify", "[", 
      RowBox[{
       RowBox[{"(", 
        RowBox[{"E\[Pi]eint\[Eta]2", "/.", 
         RowBox[{"\[Eta]2", "\[Rule]", "\[Eta]2M"}]}], ")"}], "-", 
       RowBox[{"(", 
        RowBox[{"E\[Pi]eint\[Eta]2", "/.", 
         RowBox[{"\[Eta]2", "\[Rule]", "\[Eta]2m"}]}], ")"}]}], "]"}]}], 
    ";"}]}]}]], "Input",
 CellChangeTimes->{{3.8054700091150303`*^9, 3.8054701918675756`*^9}, {
   3.8054702258364553`*^9, 3.8054702288151937`*^9}, {3.805470260126073*^9, 
   3.8054702647742195`*^9}, {3.80547029920999*^9, 3.8054703316051874`*^9}, 
   3.8054703639308367`*^9, {3.805472535748184*^9, 3.805472565820175*^9}, {
   3.8054727061315794`*^9, 3.8054727168143005`*^9}, {3.8054762957688346`*^9, 
   3.805476296623796*^9}, {3.805476347657586*^9, 3.8054763681093435`*^9}, {
   3.805476447263705*^9, 3.805476457330902*^9}, 3.8124749464845963`*^9, {
   3.8124749856097097`*^9, 3.812475173803663*^9}, {3.8124752636890583`*^9, 
   3.812475300288745*^9}, {3.8124753311486187`*^9, 3.8124753481605487`*^9}, {
   3.8124756450846305`*^9, 3.8124756501880727`*^9}, {3.812475840058958*^9, 
   3.812475870309523*^9}, 3.812476123080628*^9},
 CellLabel->"In[3]:=",ExpressionUUID->"cdf1330f-e9ba-45f1-8a3e-1254eb87e3c4"],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{
  RowBox[{"(*", 
   RowBox[{
    RowBox[{
    "Develop", " ", "quadratic", " ", "approimxations", " ", "around", " ", 
     RowBox[{"ehat", ".", " ", "This"}], " ", "follows", " ", 
     RowBox[{"Weitzman", "'"}], "s", " ", "approach"}], ",", " ", 
    RowBox[{
    "except", " ", "that", " ", "we", " ", "allow", " ", "the", " ", "slope", 
     " ", "of", " ", "marginal", " ", "profits", " ", "to", " ", "be", " ", 
     RowBox[{"uncertain", ".", " ", "Use"}], " ", "this", " ", "model", " ", 
     "to", " ", "Compute", " ", "the", " ", "TRUE", " ", "Optimal", " ", 
     "Tax", " ", "and", " ", "associated", " ", "TRUE", " ", "measure", " ", 
     "of", " ", "ESNB"}]}], "*)"}], "\[IndentingNewLine]", 
  RowBox[{"(*", 
   RowBox[{
    RowBox[{
    "Approximate", " ", "profits", " ", "around", " ", "the", " ", "optimal", 
     " ", "standard", " ", "ehat"}], ",", " ", 
    RowBox[{
    "leaving", " ", "ehat", " ", "as", " ", "an", " ", "unspecified", " ", 
     "parameter", " ", "for", " ", 
     RowBox[{"now", ".", " ", "Note"}], " ", "that", " ", "we", " ", "do", 
     " ", "not", " ", "need", " ", "to", " ", "do", " ", "the", " ", "same", 
     " ", "for", " ", "damages", " ", "since", " ", "damages", " ", "are", 
     " ", "already", " ", "quadratic", " ", "and", " ", "are", " ", "assumed",
      " ", "to", " ", "be", " ", "deterministic"}]}], "*)"}], 
  "\[IndentingNewLine]", 
  RowBox[{
   RowBox[{
    RowBox[{"\[Pi]equad", "=", 
     RowBox[{"\[Pi]0", "+", 
      RowBox[{"\[Pi]1", "*", 
       RowBox[{"(", 
        RowBox[{"e", "-", "ehat"}], ")"}]}], "+", 
      RowBox[{"0.5", "*", "\[Pi]2", "*", 
       RowBox[{
        RowBox[{"(", 
         RowBox[{"e", "-", "ehat"}], ")"}], "^", "2"}]}]}]}], ";"}], 
   "\[IndentingNewLine]", "\[IndentingNewLine]", 
   RowBox[{"(*", 
    RowBox[{
     RowBox[{
     "specify", " ", "\[Pi]0", " ", "as", " ", "\[Pi]e", " ", "evaluated", 
      " ", "when", " ", "e"}], " ", "=", " ", "ehat"}], "*)"}], 
   "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"\[Pi]0", "=", 
     RowBox[{"\[Pi]0", "/.", 
      RowBox[{"Simplify", "[", 
       RowBox[{
        RowBox[{"Solve", "[", 
         RowBox[{
          RowBox[{
           RowBox[{"0", "\[Equal]", 
            RowBox[{"(", 
             RowBox[{"\[Pi]e", "-", "\[Pi]equad"}], ")"}]}], "/.", 
           RowBox[{"e", "\[Rule]", "ehat"}]}], ",", "\[Pi]0"}], "]"}], "[", 
        RowBox[{"[", "1", "]"}], "]"}], "]"}]}]}], ";"}], 
   "\[IndentingNewLine]", 
   RowBox[{"(*", 
    RowBox[{
     RowBox[{"Since", " ", "ehat", " ", "is", " ", "not", " ", "random"}], 
     ",", " ", 
     RowBox[{
     "we", " ", "can", " ", "take", " ", "the", " ", "expectation", " ", "of",
       " ", "\[Pi]0", " ", "by", " ", "simply", " ", "setting", " ", "the", 
      " ", "\[Eta]i", " ", "terms", " ", "to", " ", "their", " ", "expected", 
      " ", 
      RowBox[{"values", ".", " ", "This"}], " ", "expectation", " ", "is", 
      " ", "performed", " ", "in", " ", "the", " ", "next", " ", "line", " ", 
      "of", " ", "code", " ", "and", " ", "will", " ", "be", " ", "used", " ",
       "in", " ", 
      RowBox[{"Weitzman", "'"}], "s", " ", "quadratic", " ", "approximation", 
      " ", "below"}]}], "*)"}], "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"E\[Pi]0", "=", 
     RowBox[{
      RowBox[{"\[Pi]0", "/.", 
       RowBox[{"\[Eta]1", "\[Rule]", "0"}]}], "/.", 
      RowBox[{"\[Eta]2", "\[Rule]", "1"}]}]}], ";"}], "\[IndentingNewLine]", 
   "\[IndentingNewLine]", 
   RowBox[{"(*", 
    RowBox[{
     RowBox[{
     "specify", " ", "\[Pi]1", " ", "as", " ", "the", " ", "derivative", " ", 
      "of", " ", "\[Pi]e", " ", "with", " ", "respect", " ", "to", " ", "e"}],
      ",", " ", 
     RowBox[{
      RowBox[{
      "given", " ", "that", " ", "e", " ", "is", " ", "evaluated", " ", "at", 
       " ", "e"}], " ", "=", " ", "ehat"}]}], "*)"}], "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"\[Pi]1", "=", 
     RowBox[{"\[Pi]1", "/.", 
      RowBox[{"Simplify", "[", 
       RowBox[{
        RowBox[{"Solve", "[", 
         RowBox[{
          RowBox[{
           RowBox[{"0", "\[Equal]", 
            RowBox[{"D", "[", 
             RowBox[{
              RowBox[{"(", 
               RowBox[{"\[Pi]e", "-", "\[Pi]equad"}], ")"}], ",", "e"}], 
             "]"}]}], "/.", 
           RowBox[{"e", "\[Rule]", "ehat"}]}], ",", "\[Pi]1"}], "]"}], "[", 
        RowBox[{"[", "1", "]"}], "]"}], "]"}]}]}], ";"}], 
   "\[IndentingNewLine]", 
   RowBox[{"(*", 
    RowBox[{
    "take", " ", "the", " ", "expectation", " ", "of", " ", "\[Pi]1", " ", 
     "to", " ", "be", " ", "used", " ", "in", " ", 
     RowBox[{"Weitzman", "'"}], "s", " ", "quadratic", " ", "approximation", 
     " ", "below"}], "*)"}], "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"E\[Pi]1", "=", 
     RowBox[{
      RowBox[{"\[Pi]1", "/.", 
       RowBox[{"\[Eta]1", "\[Rule]", "0"}]}], "/.", 
      RowBox[{"\[Eta]2", "\[Rule]", "1"}]}]}], ";"}], "\[IndentingNewLine]", 
   RowBox[{"(*", 
    RowBox[{
     RowBox[{"specify", " ", 
      RowBox[{"Weitzman", "'"}], "s", " ", "\[Alpha]", " ", "term", " ", "as",
       " ", "the", " ", "difference", " ", "between", " ", "\[Pi]1", " ", 
      "and", " ", "its", " ", "expectation"}], ",", " ", 
     RowBox[{
      RowBox[{
      "thereby", " ", "ensuring", " ", "the", " ", "coefficient", " ", "of", 
       " ", "e", " ", "is", " ", "\[Pi]1", " ", "and", " ", "that", " ", "E", 
       RowBox[{"{", "\[Alpha]", "}"}]}], " ", "=", " ", "0"}]}], "*)"}], 
   "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"\[Alpha]", "=", 
     RowBox[{"\[Pi]1", "-", "E\[Pi]1"}]}], ";"}], "\[IndentingNewLine]", 
   "\[IndentingNewLine]", 
   RowBox[{"(*", 
    RowBox[{
     RowBox[{
     "specify", " ", "\[Pi]2", " ", "as", " ", "the", " ", "second", " ", 
      "derivative", " ", "of", " ", "\[Pi]e", " ", "with", " ", "respect", 
      " ", "to", " ", "e"}], ",", " ", 
     RowBox[{
      RowBox[{
      "given", " ", "that", " ", "e", " ", "is", " ", "evaluated", " ", "at", 
       " ", "e"}], " ", "=", " ", "ehat"}]}], "*)"}], "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"\[Pi]2", "=", 
     RowBox[{"\[Pi]2", "/.", 
      RowBox[{"Simplify", "[", 
       RowBox[{
        RowBox[{"Solve", "[", 
         RowBox[{
          RowBox[{
           RowBox[{"0", "\[Equal]", 
            RowBox[{"D", "[", 
             RowBox[{
              RowBox[{"D", "[", 
               RowBox[{
                RowBox[{"(", 
                 RowBox[{"\[Pi]e", "-", "\[Pi]equad"}], ")"}], ",", "e"}], 
               "]"}], ",", "e"}], "]"}]}], "/.", 
           RowBox[{"e", "\[Rule]", "ehat"}]}], ",", "\[Pi]2"}], "]"}], "[", 
        RowBox[{"[", "1", "]"}], "]"}], "]"}]}]}], ";"}], 
   "\[IndentingNewLine]", 
   RowBox[{"(*", 
    RowBox[{
    "take", " ", "the", " ", "expectation", " ", "of", " ", "\[Pi]2", " ", 
     "to", " ", "be", " ", "used", " ", "in", " ", 
     RowBox[{"Weitzman", "'"}], "s", " ", "quadratic", " ", "approximation", 
     " ", "below"}], "*)"}], "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"E\[Pi]2", "=", 
     RowBox[{"\[Pi]2", "/.", 
      RowBox[{"\[Eta]2", "\[Rule]", " ", "1"}]}]}], ";"}], 
   "\[IndentingNewLine]", 
   RowBox[{"(*", 
    RowBox[{
     RowBox[{"specify", " ", 
      RowBox[{"Weitzman", "'"}], "s", " ", "\[Gamma]", " ", "term", " ", "as",
       " ", "the", " ", "ratio", " ", "of", " ", "\[Pi]2", " ", "to", " ", 
      "its", " ", "expectation"}], ",", " ", 
     RowBox[{
      RowBox[{
      "thereby", " ", "ensuring", " ", "the", " ", "coefficient", " ", "of", 
       " ", 
       RowBox[{"e", "^", "2"}], " ", "is", " ", "\[Pi]2", " ", "that", " ", 
       "E", 
       RowBox[{"{", "\[Gamma]", "}"}]}], " ", "=", " ", "1"}]}], "*)"}], 
   "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"\[Gamma]", "=", 
     RowBox[{"\[Pi]2", "/", "E\[Pi]2"}]}], ";"}], "\[IndentingNewLine]", 
   "\[IndentingNewLine]", 
   RowBox[{"(*", 
    RowBox[{
     RowBox[{
     "The", " ", "TRUE", " ", "quadratic", " ", "approximation", " ", "of", 
      " ", "restricted", " ", "profits"}], ",", " ", 
     RowBox[{
     "given", " ", "that", " ", "\[Gamma]", " ", "is", " ", "uncertain", " ", 
      "and", " ", "correlated", " ", "with", " ", "\[Alpha]"}]}], "*)"}], 
   "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"\[Pi]quad", "=", 
     RowBox[{"E\[Pi]0", "+", 
      RowBox[{
       RowBox[{"(", 
        RowBox[{"E\[Pi]1", "+", "\[Alpha]"}], ")"}], "*", 
       RowBox[{"(", 
        RowBox[{"e", "-", "ehat"}], ")"}]}], "+", 
      RowBox[{"0.5", "*", "E\[Pi]2", "*", "\[Gamma]", "*", 
       RowBox[{
        RowBox[{"(", 
         RowBox[{"e", "-", "ehat"}], ")"}], "^", "2"}]}]}]}], ";"}], 
   "\[IndentingNewLine]", "\[IndentingNewLine]", 
   RowBox[{"(*", 
    RowBox[{
     RowBox[{
     "Covariance", " ", "between", " ", "\[Alpha]", " ", "and", " ", 
      "\[Gamma]"}], ",", " ", 
     RowBox[{
     "as", " ", "derived", " ", "in", " ", "equation", " ", "15", " ", "in", 
      " ", "the", " ", "text"}]}], "*)"}], "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"cov\[Alpha]\[Gamma]", "=", 
     RowBox[{"E\[Pi]2", "*", "ehat", "*", "\[Sigma]2\[Eta]2"}]}], ";"}], 
   "\[IndentingNewLine]", 
   RowBox[{"(*", 
    RowBox[{
    "variance", " ", "of", " ", "\[Eta]2", " ", "given", " ", "a", " ", 
     "uniform", " ", "distribution"}], "*)"}], "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"\[Sigma]2\[Eta]2", "=", 
     RowBox[{
      RowBox[{"(", 
       RowBox[{
        RowBox[{"(", 
         RowBox[{"\[Eta]2M", "-", "\[Eta]2m"}], ")"}], "^", "2"}], " ", ")"}],
       "/", "12"}]}], ";"}], "\[IndentingNewLine]", "\[IndentingNewLine]", 
   "\[IndentingNewLine]", 
   RowBox[{"(*", 
    RowBox[{
     RowBox[{"Compute", " ", "the", " ", 
      RowBox[{"firm", "'"}], "s", " ", "response", " ", "to", " ", "the", " ",
       "tax"}], ",", " ", 
     RowBox[{"i", ".", "e", "."}], ",", " ", 
     RowBox[{"its", " ", "response", " ", "function"}], ",", " ", 
     RowBox[{
     "based", " ", "on", " ", "the", " ", "quadratic", " ", "approximation", 
      " ", "around", " ", "ehat"}]}], "*)"}], "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"e\[Tau]", "=", 
     RowBox[{"e", "/.", 
      RowBox[{
       RowBox[{"Solve", "[", 
        RowBox[{
         RowBox[{
          RowBox[{"D", "[", 
           RowBox[{
            RowBox[{"\[Pi]quad", "-", 
             RowBox[{"\[Tau]", "*", "e"}]}], ",", "e"}], "]"}], "\[Equal]", 
          "0"}], ",", "e"}], "]"}], "[", 
       RowBox[{"[", "1", "]"}], "]"}]}]}], ";"}], "\[IndentingNewLine]", 
   "\[IndentingNewLine]", 
   RowBox[{"(*", 
    RowBox[{
     RowBox[{"Define", " ", "the", " ", "optimal"}], ",", " ", 
     RowBox[{
     "approximated", " ", "profit", " ", "function", " ", "given", " ", "the",
       " ", 
      RowBox[{"firm", "'"}], "s", " ", "response", " ", "function"}]}], 
    "*)"}], " ", "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"\[Pi]\[Tau]", "=", 
     RowBox[{"\[Pi]quad", "/.", 
      RowBox[{"e", "\[Rule]", "e\[Tau]"}]}]}], ";"}], "\[IndentingNewLine]", 
   "\[IndentingNewLine]", 
   RowBox[{"(*", 
    RowBox[{
    "Caculate", " ", "expected", " ", "value", " ", "of", " ", "approximated",
      " ", "profits", " ", "given", " ", "the", " ", 
     RowBox[{"firm", "'"}], "s", " ", "response", " ", "to", " ", "the", " ", 
     RowBox[{"tax", ".", " ", "This"}], " ", "is", " ", "done", " ", "in", 
     " ", "the", " ", "same", " ", "manner", " ", "as", " ", "above"}], 
    "*)"}], " ", "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"E\[Pi]\[Tau]int\[Eta]1", "=", 
     RowBox[{
      RowBox[{"Integrate", "[", 
       RowBox[{"\[Pi]\[Tau]", ",", "\[Eta]1"}], "]"}], "/", 
      RowBox[{"(", 
       RowBox[{"\[Eta]1M", "-", "\[Eta]1m"}], ")"}]}]}], ";"}], 
   "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"E\[Pi]\[Tau]\[Eta]1", "=", 
     RowBox[{"Simplify", "[", 
      RowBox[{
       RowBox[{"(", 
        RowBox[{"E\[Pi]\[Tau]int\[Eta]1", "/.", 
         RowBox[{"\[Eta]1", "\[Rule]", "\[Eta]1M"}]}], ")"}], "-", 
       RowBox[{"(", 
        RowBox[{"E\[Pi]\[Tau]int\[Eta]1", "/.", 
         RowBox[{"\[Eta]1", "\[Rule]", "\[Eta]1m"}]}], ")"}]}], "]"}]}], 
    ";"}], "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"E\[Pi]\[Tau]int\[Eta]2", "=", 
     RowBox[{
      RowBox[{"Integrate", "[", 
       RowBox[{"E\[Pi]\[Tau]\[Eta]1", ",", "\[Eta]2"}], "]"}], "/", 
      RowBox[{"(", 
       RowBox[{"\[Eta]2M", "-", "\[Eta]2m"}], ")"}]}]}], ";"}], 
   "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"E\[Pi]\[Tau]", "=", 
     RowBox[{"Simplify", "[", 
      RowBox[{
       RowBox[{"(", 
        RowBox[{"E\[Pi]\[Tau]int\[Eta]2", "/.", 
         RowBox[{"\[Eta]2", "\[Rule]", "\[Eta]2M"}]}], ")"}], "-", 
       RowBox[{"(", 
        RowBox[{"E\[Pi]\[Tau]int\[Eta]2", "/.", 
         RowBox[{"\[Eta]2", "\[Rule]", "\[Eta]2m"}]}], ")"}]}], "]"}]}], 
    ";"}], "\[IndentingNewLine]", "\[IndentingNewLine]", 
   RowBox[{"(*", 
    RowBox[{
    "Caculate", " ", "expected", " ", "value", " ", "of", " ", "damages", " ",
      "given", " ", "the", " ", 
     RowBox[{"firm", "'"}], "s", " ", "response", " ", "to", " ", "the", " ", 
     RowBox[{"tax", ".", " ", "This"}], " ", "is", " ", "done", " ", "in", 
     " ", "the", " ", "same", " ", "manner", " ", "as", " ", "above"}], 
    "*)"}], " ", "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"dam\[Tau]", "=", 
     RowBox[{"dam", "/.", 
      RowBox[{"e", "\[Rule]", "e\[Tau]"}]}]}], ";"}], "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"Edam\[Tau]int\[Eta]1", "=", 
     RowBox[{
      RowBox[{"Integrate", "[", 
       RowBox[{"dam\[Tau]", ",", "\[Eta]1"}], "]"}], "/", 
      RowBox[{"(", 
       RowBox[{"\[Eta]1M", "-", "\[Eta]1m"}], ")"}]}]}], ";"}], 
   "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"Edam\[Tau]\[Eta]1", "=", 
     RowBox[{"Simplify", "[", 
      RowBox[{
       RowBox[{"(", 
        RowBox[{"Edam\[Tau]int\[Eta]1", "/.", 
         RowBox[{"\[Eta]1", "\[Rule]", "\[Eta]1M"}]}], ")"}], "-", 
       RowBox[{"(", 
        RowBox[{"Edam\[Tau]int\[Eta]1", "/.", 
         RowBox[{"\[Eta]1", "\[Rule]", "\[Eta]1m"}]}], ")"}]}], "]"}]}], 
    ";"}], "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"Edam\[Tau]int\[Eta]2", "=", 
     RowBox[{
      RowBox[{"Integrate", "[", 
       RowBox[{"Edam\[Tau]\[Eta]1", ",", "\[Eta]2"}], "]"}], "/", 
      RowBox[{"(", 
       RowBox[{"\[Eta]2M", "-", "\[Eta]2m"}], ")"}]}]}], ";"}], 
   "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"Edam\[Tau]", "=", 
     RowBox[{"Simplify", "[", 
      RowBox[{
       RowBox[{"(", 
        RowBox[{"Edam\[Tau]int\[Eta]2", "/.", 
         RowBox[{"\[Eta]2", "\[Rule]", "\[Eta]2M"}]}], ")"}], "-", 
       RowBox[{"(", 
        RowBox[{"Edam\[Tau]int\[Eta]2", "/.", 
         RowBox[{"\[Eta]2", "\[Rule]", "\[Eta]2m"}]}], ")"}]}], "]"}]}], 
    ";"}], "\[IndentingNewLine]", "\[IndentingNewLine]", 
   RowBox[{"(*", 
    RowBox[{
     RowBox[{
     "Solve", " ", "for", " ", "the", " ", "Truly", " ", "Optimal", " ", 
      "Tax"}], ",", " ", 
     RowBox[{"given", " ", 
      RowBox[{"ehat", ".", " ", "Mathematica"}], " ", "will", " ", "give", 
      " ", "a", " ", "warning", " ", "that", " ", "it", " ", "numericized", 
      " ", "the", " ", "result"}], ",", " ", 
     RowBox[{
     "but", " ", "the", " ", "numerical", " ", "result", " ", "will", " ", 
      "be", " ", "the", " ", "same", " ", "as", " ", "if", " ", "we", " ", 
      "were", " ", "to", " ", "have", " ", "already", " ", "specified", " ", 
      "a", " ", "value", " ", "for", " ", "ehat"}]}], "*)"}], 
   "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"\[Tau]star", "=", 
     RowBox[{"\[Tau]", "/.", 
      RowBox[{
       RowBox[{"Solve", "[", 
        RowBox[{
         RowBox[{
          RowBox[{"D", "[", 
           RowBox[{
            RowBox[{"E\[Pi]\[Tau]", "-", "Edam\[Tau]"}], ",", "\[Tau]"}], 
           "]"}], "\[Equal]", "0"}], ",", "\[Tau]"}], "]"}], "[", 
       RowBox[{"[", "1", "]"}], "]"}]}]}], ";"}], 
   "\[IndentingNewLine]"}]}]], "Input",
 CellChangeTimes->CompressedData["
1:eJxTTMoPSmViYGAQBmIQ7b78p0VM1BvHpv+m1iD6lnniWRD9L9joHIiWebba
JRZIX1U0cAXREvda7Hakv3F8WvPfAURPbT8fCqIv3X8RB6JNHqUmg+h3BWfT
QXRf7oECEK1zMbYMRIc43akB0TXeVyaC6AKj43NBtJNZ7T4QvSqq5RiIfrLz
5SkQXfd7xzkQvdeX7TzYPK+AK2B9L7ZeA9Eqj0UfgunLci9BdHe9Xt9OIH1y
4c6JILpES812F5B+E94Ppv0m73AC0TP2/XcD0d+a33iD6BNrX/uA6CTWn7bv
gPSzMEc7EC1m/icQRO9dcjYIRAMAfZipaA==
  "],
 CellLabel->"In[7]:=",ExpressionUUID->"3a775e5e-012d-4f93-9d76-951eebbda49e"],

Cell[BoxData[
 TemplateBox[{
  "Solve","ratnz",
   "\"Solve was unable to solve the system with inexact coefficients. The \
answer was obtained by solving a corresponding exact system and numericizing \
the result.\"",2,30,1,28351534326494024203,"Local"},
  "MessageTemplate"]], "Message", "MSG",
 CellChangeTimes->{3.8124798973504543`*^9, 3.812707991782346*^9, 
  3.812712046517466*^9, 3.8127121304927197`*^9, 3.8127135203579946`*^9},
 CellLabel->
  "During evaluation of \
In[7]:=",ExpressionUUID->"ff7be5fb-cd86-40ff-829b-8fa8e66dc96e"]
}, Open  ]],

Cell[BoxData[
 RowBox[{
  RowBox[{"(*", 
   RowBox[{
    RowBox[{"Construct", " ", 
     RowBox[{"Weitzman", "'"}], "s", " ", "analysis", " ", "based", " ", "on",
      " ", "a", " ", "quadratic", " ", "approximation", " ", "of", " ", "the",
      " ", "profit", " ", "function"}], ",", " ", 
    RowBox[{
     RowBox[{"given", " ", "ehat", " ", "and", " ", "with", " ", "\[Gamma]"}],
      "=", "1"}]}], "*)"}], "\n", 
  RowBox[{
   RowBox[{
    RowBox[{"\[Pi]W", "=", 
     RowBox[{"E\[Pi]0", "+", 
      RowBox[{
       RowBox[{"(", 
        RowBox[{"E\[Pi]1", "+", "\[Alpha]"}], ")"}], "*", 
       RowBox[{"(", 
        RowBox[{"e", "-", "ehat"}], ")"}]}], "+", 
      RowBox[{"0.5", "*", "E\[Pi]2", "*", 
       RowBox[{
        RowBox[{"(", 
         RowBox[{"e", "-", "ehat"}], ")"}], "^", "2"}]}]}]}], ";"}], 
   "\[IndentingNewLine]", "\[IndentingNewLine]", 
   RowBox[{"(*", 
    RowBox[{
     RowBox[{"The", " ", 
      RowBox[{"firm", "'"}], "s", " ", "response", " ", "to", " ", "the", " ",
       "tax"}], ",", " ", 
     RowBox[{
      RowBox[{
      "given", " ", "ehat", " ", "and", " ", "that", " ", "\[Gamma]"}], "=", 
      "1"}]}], "*)"}], "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"e\[Tau]W", "=", 
     RowBox[{"e", "/.", 
      RowBox[{
       RowBox[{"Solve", "[", 
        RowBox[{
         RowBox[{"D", "[", 
          RowBox[{
           RowBox[{
            RowBox[{"\[Pi]W", "-", 
             RowBox[{"\[Tau]W", "*", "e"}]}], "\[Equal]", "0"}], ",", "e"}], 
          "]"}], ",", "e"}], "]"}], "[", 
       RowBox[{"[", "1", "]"}], "]"}]}]}], ";"}], "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"\[Pi]\[Tau]W", "=", 
     RowBox[{"\[Pi]W", "/.", 
      RowBox[{"e", "\[Rule]", "e\[Tau]W"}]}]}], ";"}], "\[IndentingNewLine]", 
   "\[IndentingNewLine]", 
   RowBox[{"(*", 
    RowBox[{
    "Caculate", " ", "expected", " ", "value", " ", "of", " ", "approximated",
      " ", "profits", " ", "given", " ", "the", " ", 
     RowBox[{"firm", "'"}], "s", " ", "response", " ", "to", " ", "the", " ", 
     RowBox[{"tax", ".", " ", "This"}], " ", "is", " ", "done", " ", "in", 
     " ", "the", " ", "same", " ", "manner", " ", "as", " ", "above"}], 
    "*)"}], " ", "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"E\[Pi]\[Tau]Wint\[Eta]1", "=", 
     RowBox[{
      RowBox[{"Integrate", "[", 
       RowBox[{"\[Pi]\[Tau]W", ",", "\[Eta]1"}], "]"}], "/", 
      RowBox[{"(", 
       RowBox[{"\[Eta]1M", "-", "\[Eta]1m"}], ")"}]}]}], ";"}], 
   "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"E\[Pi]\[Tau]W\[Eta]1", "=", 
     RowBox[{"Simplify", "[", 
      RowBox[{
       RowBox[{"(", 
        RowBox[{"E\[Pi]\[Tau]Wint\[Eta]1", "/.", 
         RowBox[{"\[Eta]1", "\[Rule]", "\[Eta]1M"}]}], ")"}], "-", 
       RowBox[{"(", 
        RowBox[{"E\[Pi]\[Tau]Wint\[Eta]1", "/.", 
         RowBox[{"\[Eta]1", "\[Rule]", "\[Eta]1m"}]}], ")"}]}], "]"}]}], 
    ";"}], "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"E\[Pi]\[Tau]Wint\[Eta]2", "=", 
     RowBox[{
      RowBox[{"Integrate", "[", 
       RowBox[{"E\[Pi]\[Tau]W\[Eta]1", ",", "\[Eta]2"}], "]"}], "/", 
      RowBox[{"(", 
       RowBox[{"\[Eta]2M", "-", "\[Eta]2m"}], ")"}]}]}], ";"}], 
   "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"E\[Pi]\[Tau]W", "=", 
     RowBox[{"Simplify", "[", 
      RowBox[{
       RowBox[{"(", 
        RowBox[{"E\[Pi]\[Tau]Wint\[Eta]2", "/.", 
         RowBox[{"\[Eta]2", "\[Rule]", "\[Eta]2M"}]}], ")"}], "-", 
       RowBox[{"(", 
        RowBox[{"E\[Pi]\[Tau]Wint\[Eta]2", "/.", 
         RowBox[{"\[Eta]2", "\[Rule]", "\[Eta]2m"}]}], ")"}]}], "]"}]}], 
    ";"}], "\[IndentingNewLine]", "\[IndentingNewLine]", 
   RowBox[{"(*", 
    RowBox[{
    "Caculate", " ", "expected", " ", "value", " ", "of", " ", "damages", " ",
      "given", " ", "the", " ", 
     RowBox[{"firm", "'"}], "s", " ", "response", " ", "to", " ", "the", " ", 
     RowBox[{"tax", ".", " ", "This"}], " ", "is", " ", "done", " ", "in", 
     " ", "the", " ", "same", " ", "manner", " ", "as", " ", "above"}], 
    "*)"}], " ", "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"dam\[Tau]W", "=", 
     RowBox[{"dam", "/.", 
      RowBox[{"e", "\[Rule]", "e\[Tau]W"}]}]}], ";"}], "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"Edam\[Tau]Wint\[Eta]1", "=", 
     RowBox[{
      RowBox[{"Integrate", "[", 
       RowBox[{"dam\[Tau]W", ",", "\[Eta]1"}], "]"}], "/", 
      RowBox[{"(", 
       RowBox[{"\[Eta]1M", "-", "\[Eta]1m"}], ")"}]}]}], ";"}], 
   "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"Edam\[Tau]W\[Eta]1", "=", 
     RowBox[{"Simplify", "[", 
      RowBox[{
       RowBox[{"(", 
        RowBox[{"Edam\[Tau]Wint\[Eta]1", "/.", 
         RowBox[{"\[Eta]1", "\[Rule]", "\[Eta]1M"}]}], ")"}], "-", 
       RowBox[{"(", 
        RowBox[{"Edam\[Tau]Wint\[Eta]1", "/.", 
         RowBox[{"\[Eta]1", "\[Rule]", "\[Eta]1m"}]}], ")"}]}], "]"}]}], 
    ";"}], "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"Edam\[Tau]Wint\[Eta]2", "=", 
     RowBox[{
      RowBox[{"Integrate", "[", 
       RowBox[{"Edam\[Tau]W\[Eta]1", ",", "\[Eta]2"}], "]"}], "/", 
      RowBox[{"(", 
       RowBox[{"\[Eta]2M", "-", "\[Eta]2m"}], ")"}]}]}], ";"}], 
   "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"Edam\[Tau]W", "=", 
     RowBox[{"Simplify", "[", 
      RowBox[{
       RowBox[{"(", 
        RowBox[{"Edam\[Tau]Wint\[Eta]2", "/.", 
         RowBox[{"\[Eta]2", "\[Rule]", "\[Eta]2M"}]}], ")"}], "-", 
       RowBox[{"(", 
        RowBox[{"Edam\[Tau]Wint\[Eta]2", "/.", 
         RowBox[{"\[Eta]2", "\[Rule]", "\[Eta]2m"}]}], ")"}]}], "]"}]}], 
    ";"}], "\[IndentingNewLine]"}]}]], "Input",
 CellChangeTimes->{{3.8054781857205915`*^9, 3.805478201229881*^9}, {
   3.8054782323862457`*^9, 3.805478255877183*^9}, {3.805478294839895*^9, 
   3.8054784702894406`*^9}, {3.805478514863439*^9, 3.80547855238476*^9}, {
   3.8054785840141525`*^9, 3.8054785935257053`*^9}, 3.805480406985426*^9, {
   3.812479600224328*^9, 3.8124797221226287`*^9}, {3.8127117481239796`*^9, 
   3.812711819325647*^9}, 3.8127121211908903`*^9, 3.8127121861552424`*^9},
 CellLabel->"In[31]:=",ExpressionUUID->"701aced9-601e-49c4-9f34-1b9fcca91619"],

Cell[BoxData[
 RowBox[{
  RowBox[{"(*", 
   RowBox[{
   "Compute", " ", "the", " ", "Optimal", " ", "Standard", " ", "and", " ", 
    "marginal", " ", "damages", " ", "using", " ", "the", " ", "original", 
    " ", 
    RowBox[{"model", ".", " ", "The"}], " ", "solution", " ", "will", " ", 
    "be", " ", "an", " ", "analytical", " ", "one", " ", "until", " ", "we", 
    " ", "specify", " ", "parameter", " ", "values"}], "*)"}], 
  "\[IndentingNewLine]", 
  RowBox[{
   RowBox[{
    RowBox[{"ehat", "=", 
     RowBox[{"e", "/.", 
      RowBox[{
       RowBox[{"Simplify", "[", 
        RowBox[{"Solve", "[", 
         RowBox[{
          RowBox[{
           RowBox[{"D", "[", 
            RowBox[{
             RowBox[{"E\[Pi]e", "-", "dam"}], ",", "e"}], "]"}], "\[Equal]", 
           "0"}], ",", "e"}], "]"}], "]"}], "[", 
       RowBox[{"[", "1", "]"}], "]"}]}]}], ";"}], "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"MD", " ", "=", " ", 
     RowBox[{
      RowBox[{"D", "[", 
       RowBox[{"dam", ",", "e"}], "]"}], "/.", 
      RowBox[{"e", "\[Rule]", "ehat"}]}]}], ";"}]}]}]], "Input",
 CellChangeTimes->{{3.8054704252426195`*^9, 3.8054704663591633`*^9}, {
   3.805470511678501*^9, 3.805470650102188*^9}, 3.805474674597808*^9, {
   3.8054752351169806`*^9, 3.8054752400709543`*^9}, {3.8054776871272364`*^9, 
   3.805477687524078*^9}, {3.8054794997891607`*^9, 3.805479500590164*^9}, {
   3.8124797437138457`*^9, 3.812479805175641*^9}, {3.812710275458621*^9, 
   3.812710329754256*^9}},
 CellLabel->"In[43]:=",ExpressionUUID->"eded3f86-e60f-4f25-a0af-902672acae05"],

Cell[BoxData[
 RowBox[{
  RowBox[{"(*", 
   RowBox[{"Parameter", " ", "values"}], "*)"}], "\[IndentingNewLine]", 
  RowBox[{
   RowBox[{
    RowBox[{"\[Pi]e0", "=", "0"}], ";"}], "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"\[Pi]e1", "=", "1000"}], ";"}], "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"\[Pi]e2", "=", 
     RowBox[{"-", "100"}]}], ";"}], "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"d1", "=", "45"}], ";"}], "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"\[Eta]1m", "=", 
     RowBox[{
      RowBox[{"-", "0.1"}], "*", "\[Pi]e1"}]}], ";"}], "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"\[Eta]1M", "=", 
     RowBox[{"-", "\[Eta]1m"}]}], ";"}], "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"\[Eta]2m", "=", "0.1"}], ";"}], "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"\[Eta]2M", "=", 
     RowBox[{"2", "-", " ", "\[Eta]2m"}]}], ";"}]}]}]], "Input",
 CellLabel->"In[45]:=",ExpressionUUID->"e26b58c6-972a-48b9-8b6d-3ed929996a4b"],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"\[IndentingNewLine]", "\[IndentingNewLine]", 
  RowBox[{"(*", 
   RowBox[{"Numerical", " ", "Results"}], "*)"}], "\[IndentingNewLine]", 
  RowBox[{"\"\<Unregulated ESNB\>\"", "\[IndentingNewLine]", 
   RowBox[{"ESNB0", "=", 
    RowBox[{
     RowBox[{"(", 
      RowBox[{"E\[Pi]\[Tau]", "-", "Edam\[Tau]"}], ")"}], "/.", 
     RowBox[{"\[Tau]", "\[Rule]", "0"}]}]}], "\[IndentingNewLine]", 
   "\[IndentingNewLine]", 
   "\"\<% Reduction in emissions from baseline under optimal standard\>\"", 
   "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{
     RowBox[{"-", "ehat"}], "/", 
     RowBox[{"(", 
      RowBox[{"\[Pi]e1", "/", "\[Pi]e2"}], ")"}]}], "*", "100"}], 
   "\[IndentingNewLine]", "\[IndentingNewLine]", "\"\<ESNB e hat\>\"", 
   "\[IndentingNewLine]", 
   RowBox[{"ESNBhat", "=", 
    RowBox[{
     RowBox[{"(", 
      RowBox[{"E\[Pi]e", "-", "dam"}], ")"}], "/.", 
     RowBox[{"e", "\[Rule]", "ehat"}]}]}], " ", "\[IndentingNewLine]", 
   "\[IndentingNewLine]", 
   "\"\<Expected Gain in SNB from regulating under the standard\>\"", 
   "\[IndentingNewLine]", 
   RowBox[{"ESNBhat", "-", "ESNB0"}], "\[IndentingNewLine]", 
   "\[IndentingNewLine]", "\"\<Covariance between \[Alpha] and \[Gamma]\>\"", 
   "\[IndentingNewLine]", "cov\[Alpha]\[Gamma]", "\[IndentingNewLine]", 
   "\[IndentingNewLine]", "\"\<true tax\>\"", "\[IndentingNewLine]", 
   "\[Tau]star", "\[IndentingNewLine]", "\[IndentingNewLine]", "\"\<r\>\"", 
   "\[IndentingNewLine]", 
   RowBox[{"\[Tau]star", "-", "MD"}], "\[IndentingNewLine]", 
   "\[IndentingNewLine]", "\"\<ESNB True \[Tau]\>\"", "\[IndentingNewLine]", 
   RowBox[{"ESNB\[Tau]", "=", 
    RowBox[{
     RowBox[{"(", 
      RowBox[{"E\[Pi]\[Tau]", "-", "Edam\[Tau]"}], ")"}], "/.", 
     RowBox[{"\[Tau]", "\[Rule]", "\[Tau]star"}]}]}], "\[IndentingNewLine]", 
   "\[IndentingNewLine]", 
   "\"\<Expected Gain in SNB from regulating under the true tax\>\"", 
   "\[IndentingNewLine]", 
   RowBox[{"ESNB\[Tau]", "-", "ESNB0"}], "\[IndentingNewLine]", 
   "\[IndentingNewLine]", "\"\<True \[CapitalDelta]\>\"", 
   "\[IndentingNewLine]", 
   RowBox[{"\[CapitalDelta]opt", "=", 
    RowBox[{"ESNB\[Tau]", "-", "ESNBhat"}]}], "\[IndentingNewLine]", 
   "\[IndentingNewLine]", 
   "\"\<Expected SNB for an unregulated system, assuming \[Gamma]=1\>\"", 
   "\[IndentingNewLine]", 
   RowBox[{"ESNBW0", "=", 
    RowBox[{
     RowBox[{"(", 
      RowBox[{"E\[Pi]\[Tau]W", "-", "Edam\[Tau]W"}], ")"}], "/.", 
     RowBox[{"\[Tau]W", "\[Rule]", "0"}]}]}], "\[IndentingNewLine]", "\n", 
   "\"\<Weitzman's Tax = MD at ehat\>\"", "\[IndentingNewLine]", 
   RowBox[{"\[Tau]W", "=", 
    RowBox[{"2", "*", "d1", "*", "ehat"}]}], "\[IndentingNewLine]", 
   "\[IndentingNewLine]", "\"\<ESNB Tax Weitman\>\"", "\[IndentingNewLine]", 
   RowBox[{"ESNB\[Tau]W", " ", "=", " ", 
    RowBox[{
     RowBox[{"(", 
      RowBox[{"E\[Pi]\[Tau]W", "-", "Edam\[Tau]W"}], ")"}], "/.", 
     RowBox[{"\[Tau]", "\[Rule]", "\[Tau]W"}]}]}], "\[IndentingNewLine]", 
   "\[IndentingNewLine]", 
   "\"\<Expected Gain in SNB from regulating under Weitzman's tax and \
assuming \[Gamma]=1 when deriving unregulated ESNB\>\"", 
   "\[IndentingNewLine]", 
   RowBox[{"ESNB\[Tau]W", "-", "ESNBW0"}], "\[IndentingNewLine]", 
   "\[IndentingNewLine]", "\"\<\[CapitalDelta]W\>\"", "\[IndentingNewLine]", 
   " ", 
   RowBox[{"ESNB\[Tau]W", "-", "ESNBhat"}], "\[IndentingNewLine]", 
   "\[IndentingNewLine]", "\"\<relative slopes at ehat\>\"", 
   "\[IndentingNewLine]", 
   RowBox[{"(", 
    RowBox[{
     RowBox[{"-", "E\[Pi]2"}], "-", 
     RowBox[{"2", "*", "d1"}]}], ")"}], "\[IndentingNewLine]", 
   "\[IndentingNewLine]", "\[IndentingNewLine]", "\[IndentingNewLine]", 
   "\"\<ESNB True from using Weitzman's tax\>\"", "\[IndentingNewLine]", 
   RowBox[{"ESNB\[Tau]Wtrue", "=", 
    RowBox[{
     RowBox[{"(", 
      RowBox[{"E\[Pi]\[Tau]", "-", "Edam\[Tau]"}], ")"}], "/.", 
     RowBox[{"\[Tau]", "\[Rule]", "\[Tau]W"}]}]}], "\[IndentingNewLine]", 
   "\[IndentingNewLine]", 
   "\"\<True Expected Gain in SNB from regulating under Weitzman's tax\>\"", 
   "\[IndentingNewLine]", 
   RowBox[{"Egaintrue", "=", 
    RowBox[{"ESNB\[Tau]Wtrue", "-", "ESNB0"}]}], "\[IndentingNewLine]", 
   "\[IndentingNewLine]", "\[IndentingNewLine]", 
   "\"\<% diff in ESNB from policy based on correlated uncertainty\>\"", 
   "\[IndentingNewLine]", 
   RowBox[{"100", "*", 
    RowBox[{
     RowBox[{"(", 
      RowBox[{
       RowBox[{"Max", "[", 
        RowBox[{"ESNBhat", ",", "ESNB\[Tau]"}], "]"}], "-", 
       "ESNB\[Tau]Wtrue"}], ")"}], "/", 
     RowBox[{"Max", "[", 
      RowBox[{"ESNBhat", ",", "ESNB\[Tau]"}], "]"}]}]}], 
   "\[IndentingNewLine]", "\[IndentingNewLine]", 
   "\"\<% diff in gain in ESNB from policy based on correlated \
uncertainty\>\"", "\[IndentingNewLine]", 
   RowBox[{"100", "*", 
    RowBox[{
     RowBox[{"(", 
      RowBox[{
       RowBox[{"Max", "[", 
        RowBox[{"ESNBhat", ",", "ESNB\[Tau]"}], "]"}], "-", 
       "ESNB\[Tau]Wtrue"}], ")"}], "/", 
     RowBox[{"(", 
      RowBox[{
       RowBox[{"Max", "[", 
        RowBox[{"ESNBhat", ",", "ESNB\[Tau]"}], "]"}], "-", "ESNB0"}], 
      ")"}]}]}], "\[IndentingNewLine]", "\[IndentingNewLine]"}]}]], "Input",
 CellChangeTimes->{{3.8054776984049287`*^9, 3.8054777789768105`*^9}, {
   3.80547781960773*^9, 3.805477862935916*^9}, 3.805478812273633*^9, {
   3.80547888808515*^9, 3.805478907407302*^9}, 3.8054792744008684`*^9, {
   3.805479368060747*^9, 3.80547941494519*^9}, {3.8054797598970757`*^9, 
   3.80547976634675*^9}, 3.8054798011313214`*^9, {3.8054799193506203`*^9, 
   3.805479920276123*^9}, {3.8054799520269613`*^9, 3.80547996061756*^9}, {
   3.805480819572815*^9, 3.805480974966502*^9}, {3.805481011832762*^9, 
   3.805481138600116*^9}, {3.805481206158195*^9, 3.8054813229988976`*^9}, {
   3.805481513562931*^9, 3.805481565247733*^9}, {3.8054816120966396`*^9, 
   3.8054816601816883`*^9}, 3.8054816932465177`*^9, {3.8054823738316326`*^9, 
   3.8054825109551673`*^9}, {3.805482628993081*^9, 3.8054826975500803`*^9}, {
   3.8054829353109603`*^9, 3.8054829565867977`*^9}, {3.8054830062525845`*^9, 
   3.8054830074296217`*^9}, {3.805502375933859*^9, 3.8055023883903723`*^9}, 
   3.8055024642686963`*^9, {3.805502523630508*^9, 3.8055025613245735`*^9}, {
   3.8055025936193266`*^9, 3.8055026157269707`*^9}, {3.8055026725212507`*^9, 
   3.8055028904464884`*^9}, {3.80553909134046*^9, 3.8055390915463066`*^9}, {
   3.805539194814892*^9, 3.8055392004579077`*^9}, {3.8055392613587236`*^9, 
   3.8055392947690964`*^9}, {3.8055476373248143`*^9, 3.805547640999133*^9}, {
   3.805547698797248*^9, 3.8055476989091845`*^9}, {3.8055477429814043`*^9, 
   3.805547782204968*^9}, {3.8060612365161753`*^9, 3.806061237531761*^9}, {
   3.8060613755395637`*^9, 3.8060613980997257`*^9}, {3.806061436080139*^9, 
   3.806061484807577*^9}, {3.806061893132056*^9, 3.806061940789712*^9}, {
   3.806061987761306*^9, 3.806061987870246*^9}, {3.8060620730942707`*^9, 
   3.8060620759835773`*^9}, {3.8060621209569564`*^9, 3.806062247175007*^9}, {
   3.8124798107399144`*^9, 3.812479845425487*^9}, {3.8124805131320343`*^9, 
   3.812480551544015*^9}, {3.812708184227682*^9, 3.8127082265349627`*^9}, {
   3.812708271787597*^9, 3.812708542965426*^9}, {3.812710344409814*^9, 
   3.8127103636003323`*^9}, {3.8127104132837906`*^9, 
   3.8127104151092067`*^9}, {3.8127105446691933`*^9, 
   3.8127106025002856`*^9}, {3.812710784890956*^9, 3.8127108088113623`*^9}, {
   3.812710846254524*^9, 3.812710847475974*^9}, {3.8127108788099318`*^9, 
   3.812710901003096*^9}, {3.812710956917696*^9, 3.8127109876849403`*^9}, {
   3.8127110221249285`*^9, 3.81271105771202*^9}, {3.8127115818923197`*^9, 
   3.812711643311558*^9}, {3.812711978652289*^9, 3.8127119874960117`*^9}, {
   3.812713395376194*^9, 3.812713504567997*^9}},
 CellLabel->"In[54]:=",ExpressionUUID->"daa7a5ae-3fc1-4c4c-b9ab-118cd11e24be"],

Cell[BoxData["\<\"Unregulated ESNB\"\>"], "Output",
 CellChangeTimes->{
  3.8054797335618596`*^9, 3.805479768529723*^9, 3.805479962547964*^9, 
   3.805480011263193*^9, 3.805480158628708*^9, 3.8054805207302823`*^9, {
   3.8054806478378954`*^9, 3.8054807227787147`*^9}, {3.805480754776679*^9, 
   3.805480802785015*^9}, {3.805480860379613*^9, 3.805480975306303*^9}, {
   3.805481055037532*^9, 3.805481138991852*^9}, {3.805481207897601*^9, 
   3.8054813233177137`*^9}, {3.805481668996467*^9, 3.8054816937761745`*^9}, {
   3.8054824596940994`*^9, 3.805482511675572*^9}, 3.80548270058698*^9, {
   3.805482938776911*^9, 3.805482957449237*^9}, 3.8054830084992285`*^9, 
   3.8055023938153105`*^9, 3.805502465059903*^9, 3.8055025620124936`*^9, {
   3.8055025944080496`*^9, 3.805502624323163*^9}, {3.8055026772153587`*^9, 
   3.805502783940276*^9}, {3.805502854239046*^9, 3.80550289141193*^9}, 
   3.805539092142848*^9, 3.8055392008147535`*^9, {3.805539266359166*^9, 
   3.8055392950748816`*^9}, 3.8055476424770594`*^9, {3.805547743548753*^9, 
   3.805547782692649*^9}, 3.806061311121106*^9, {3.8060613762090282`*^9, 
   3.806061399563547*^9}, 3.806061485787599*^9, {3.8060618969884276`*^9, 
   3.8060619413353577`*^9}, 3.8060619883466463`*^9, 3.806062076635193*^9, {
   3.806062121638563*^9, 3.8060622477261877`*^9}, 3.8061436402582865`*^9, 
   3.812479908939375*^9, {3.8124805431737375`*^9, 3.8124805521521053`*^9}, 
   3.8127080260794463`*^9, 3.8127082080786448`*^9, 3.8127082767170734`*^9, 
   3.8127085542792435`*^9, 3.8127103778091154`*^9, 3.81271041590324*^9, 
   3.8127105663518934`*^9, 3.8127106029664683`*^9, 3.8127110945911274`*^9, {
   3.8127116204864073`*^9, 3.8127116438698387`*^9}, 3.8127120512537746`*^9, 
   3.8127132694632845`*^9, 3.8127135282045836`*^9},
 CellLabel->"Out[54]=",ExpressionUUID->"a027acb6-16ca-40c1-8a52-5f3b40c773be"],

Cell[BoxData[
 RowBox[{"-", "15556.897406504442`"}]], "Output",
 CellChangeTimes->{
  3.8054797335618596`*^9, 3.805479768529723*^9, 3.805479962547964*^9, 
   3.805480011263193*^9, 3.805480158628708*^9, 3.8054805207302823`*^9, {
   3.8054806478378954`*^9, 3.8054807227787147`*^9}, {3.805480754776679*^9, 
   3.805480802785015*^9}, {3.805480860379613*^9, 3.805480975306303*^9}, {
   3.805481055037532*^9, 3.805481138991852*^9}, {3.805481207897601*^9, 
   3.8054813233177137`*^9}, {3.805481668996467*^9, 3.8054816937761745`*^9}, {
   3.8054824596940994`*^9, 3.805482511675572*^9}, 3.80548270058698*^9, {
   3.805482938776911*^9, 3.805482957449237*^9}, 3.8054830084992285`*^9, 
   3.8055023938153105`*^9, 3.805502465059903*^9, 3.8055025620124936`*^9, {
   3.8055025944080496`*^9, 3.805502624323163*^9}, {3.8055026772153587`*^9, 
   3.805502783940276*^9}, {3.805502854239046*^9, 3.80550289141193*^9}, 
   3.805539092142848*^9, 3.8055392008147535`*^9, {3.805539266359166*^9, 
   3.8055392950748816`*^9}, 3.8055476424770594`*^9, {3.805547743548753*^9, 
   3.805547782692649*^9}, 3.806061311121106*^9, {3.8060613762090282`*^9, 
   3.806061399563547*^9}, 3.806061485787599*^9, {3.8060618969884276`*^9, 
   3.8060619413353577`*^9}, 3.8060619883466463`*^9, 3.806062076635193*^9, {
   3.806062121638563*^9, 3.8060622477261877`*^9}, 3.8061436402582865`*^9, 
   3.812479908939375*^9, {3.8124805431737375`*^9, 3.8124805521521053`*^9}, 
   3.8127080260794463`*^9, 3.8127082080786448`*^9, 3.8127082767170734`*^9, 
   3.8127085542792435`*^9, 3.8127103778091154`*^9, 3.81271041590324*^9, 
   3.8127105663518934`*^9, 3.8127106029664683`*^9, 3.8127110945911274`*^9, {
   3.8127116204864073`*^9, 3.8127116438698387`*^9}, 3.8127120512537746`*^9, 
   3.8127132694632845`*^9, 3.81271352820957*^9},
 CellLabel->"Out[55]=",ExpressionUUID->"9d51ae94-dd48-4adb-b153-8c131a503335"],

Cell[BoxData["\<\"% Reduction in emissions from baseline under optimal \
standard\"\>"], "Output",
 CellChangeTimes->{
  3.8054797335618596`*^9, 3.805479768529723*^9, 3.805479962547964*^9, 
   3.805480011263193*^9, 3.805480158628708*^9, 3.8054805207302823`*^9, {
   3.8054806478378954`*^9, 3.8054807227787147`*^9}, {3.805480754776679*^9, 
   3.805480802785015*^9}, {3.805480860379613*^9, 3.805480975306303*^9}, {
   3.805481055037532*^9, 3.805481138991852*^9}, {3.805481207897601*^9, 
   3.8054813233177137`*^9}, {3.805481668996467*^9, 3.8054816937761745`*^9}, {
   3.8054824596940994`*^9, 3.805482511675572*^9}, 3.80548270058698*^9, {
   3.805482938776911*^9, 3.805482957449237*^9}, 3.8054830084992285`*^9, 
   3.8055023938153105`*^9, 3.805502465059903*^9, 3.8055025620124936`*^9, {
   3.8055025944080496`*^9, 3.805502624323163*^9}, {3.8055026772153587`*^9, 
   3.805502783940276*^9}, {3.805502854239046*^9, 3.80550289141193*^9}, 
   3.805539092142848*^9, 3.8055392008147535`*^9, {3.805539266359166*^9, 
   3.8055392950748816`*^9}, 3.8055476424770594`*^9, {3.805547743548753*^9, 
   3.805547782692649*^9}, 3.806061311121106*^9, {3.8060613762090282`*^9, 
   3.806061399563547*^9}, 3.806061485787599*^9, {3.8060618969884276`*^9, 
   3.8060619413353577`*^9}, 3.8060619883466463`*^9, 3.806062076635193*^9, {
   3.806062121638563*^9, 3.8060622477261877`*^9}, 3.8061436402582865`*^9, 
   3.812479908939375*^9, {3.8124805431737375`*^9, 3.8124805521521053`*^9}, 
   3.8127080260794463`*^9, 3.8127082080786448`*^9, 3.8127082767170734`*^9, 
   3.8127085542792435`*^9, 3.8127103778091154`*^9, 3.81271041590324*^9, 
   3.8127105663518934`*^9, 3.8127106029664683`*^9, 3.8127110945911274`*^9, {
   3.8127116204864073`*^9, 3.8127116438698387`*^9}, 3.8127120512537746`*^9, 
   3.8127132694632845`*^9, 3.812713528212562*^9},
 CellLabel->"Out[56]=",ExpressionUUID->"fba2204a-19a0-47ab-aba7-be81f1c8ade0"],

Cell[BoxData["52.63157894736843`"], "Output",
 CellChangeTimes->{
  3.8054797335618596`*^9, 3.805479768529723*^9, 3.805479962547964*^9, 
   3.805480011263193*^9, 3.805480158628708*^9, 3.8054805207302823`*^9, {
   3.8054806478378954`*^9, 3.8054807227787147`*^9}, {3.805480754776679*^9, 
   3.805480802785015*^9}, {3.805480860379613*^9, 3.805480975306303*^9}, {
   3.805481055037532*^9, 3.805481138991852*^9}, {3.805481207897601*^9, 
   3.8054813233177137`*^9}, {3.805481668996467*^9, 3.8054816937761745`*^9}, {
   3.8054824596940994`*^9, 3.805482511675572*^9}, 3.80548270058698*^9, {
   3.805482938776911*^9, 3.805482957449237*^9}, 3.8054830084992285`*^9, 
   3.8055023938153105`*^9, 3.805502465059903*^9, 3.8055025620124936`*^9, {
   3.8055025944080496`*^9, 3.805502624323163*^9}, {3.8055026772153587`*^9, 
   3.805502783940276*^9}, {3.805502854239046*^9, 3.80550289141193*^9}, 
   3.805539092142848*^9, 3.8055392008147535`*^9, {3.805539266359166*^9, 
   3.8055392950748816`*^9}, 3.8055476424770594`*^9, {3.805547743548753*^9, 
   3.805547782692649*^9}, 3.806061311121106*^9, {3.8060613762090282`*^9, 
   3.806061399563547*^9}, 3.806061485787599*^9, {3.8060618969884276`*^9, 
   3.8060619413353577`*^9}, 3.8060619883466463`*^9, 3.806062076635193*^9, {
   3.806062121638563*^9, 3.8060622477261877`*^9}, 3.8061436402582865`*^9, 
   3.812479908939375*^9, {3.8124805431737375`*^9, 3.8124805521521053`*^9}, 
   3.8127080260794463`*^9, 3.8127082080786448`*^9, 3.8127082767170734`*^9, 
   3.8127085542792435`*^9, 3.8127103778091154`*^9, 3.81271041590324*^9, 
   3.8127105663518934`*^9, 3.8127106029664683`*^9, 3.8127110945911274`*^9, {
   3.8127116204864073`*^9, 3.8127116438698387`*^9}, 3.8127120512537746`*^9, 
   3.8127132694632845`*^9, 3.8127135282165537`*^9},
 CellLabel->"Out[57]=",ExpressionUUID->"f7094cb9-7693-4d26-9455-3d78a1757817"],

Cell[BoxData["\<\"ESNB e hat\"\>"], "Output",
 CellChangeTimes->{
  3.8054797335618596`*^9, 3.805479768529723*^9, 3.805479962547964*^9, 
   3.805480011263193*^9, 3.805480158628708*^9, 3.8054805207302823`*^9, {
   3.8054806478378954`*^9, 3.8054807227787147`*^9}, {3.805480754776679*^9, 
   3.805480802785015*^9}, {3.805480860379613*^9, 3.805480975306303*^9}, {
   3.805481055037532*^9, 3.805481138991852*^9}, {3.805481207897601*^9, 
   3.8054813233177137`*^9}, {3.805481668996467*^9, 3.8054816937761745`*^9}, {
   3.8054824596940994`*^9, 3.805482511675572*^9}, 3.80548270058698*^9, {
   3.805482938776911*^9, 3.805482957449237*^9}, 3.8054830084992285`*^9, 
   3.8055023938153105`*^9, 3.805502465059903*^9, 3.8055025620124936`*^9, {
   3.8055025944080496`*^9, 3.805502624323163*^9}, {3.8055026772153587`*^9, 
   3.805502783940276*^9}, {3.805502854239046*^9, 3.80550289141193*^9}, 
   3.805539092142848*^9, 3.8055392008147535`*^9, {3.805539266359166*^9, 
   3.8055392950748816`*^9}, 3.8055476424770594`*^9, {3.805547743548753*^9, 
   3.805547782692649*^9}, 3.806061311121106*^9, {3.8060613762090282`*^9, 
   3.806061399563547*^9}, 3.806061485787599*^9, {3.8060618969884276`*^9, 
   3.8060619413353577`*^9}, 3.8060619883466463`*^9, 3.806062076635193*^9, {
   3.806062121638563*^9, 3.8060622477261877`*^9}, 3.8061436402582865`*^9, 
   3.812479908939375*^9, {3.8124805431737375`*^9, 3.8124805521521053`*^9}, 
   3.8127080260794463`*^9, 3.8127082080786448`*^9, 3.8127082767170734`*^9, 
   3.8127085542792435`*^9, 3.8127103778091154`*^9, 3.81271041590324*^9, 
   3.8127105663518934`*^9, 3.8127106029664683`*^9, 3.8127110945911274`*^9, {
   3.8127116204864073`*^9, 3.8127116438698387`*^9}, 3.8127120512537746`*^9, 
   3.8127132694632845`*^9, 3.812713528218546*^9},
 CellLabel->"Out[58]=",ExpressionUUID->"62f62a72-77ff-4311-bcc5-6bec6e75c0f2"],

Cell[BoxData["2631.578947368421`"], "Output",
 CellChangeTimes->{
  3.8054797335618596`*^9, 3.805479768529723*^9, 3.805479962547964*^9, 
   3.805480011263193*^9, 3.805480158628708*^9, 3.8054805207302823`*^9, {
   3.8054806478378954`*^9, 3.8054807227787147`*^9}, {3.805480754776679*^9, 
   3.805480802785015*^9}, {3.805480860379613*^9, 3.805480975306303*^9}, {
   3.805481055037532*^9, 3.805481138991852*^9}, {3.805481207897601*^9, 
   3.8054813233177137`*^9}, {3.805481668996467*^9, 3.8054816937761745`*^9}, {
   3.8054824596940994`*^9, 3.805482511675572*^9}, 3.80548270058698*^9, {
   3.805482938776911*^9, 3.805482957449237*^9}, 3.8054830084992285`*^9, 
   3.8055023938153105`*^9, 3.805502465059903*^9, 3.8055025620124936`*^9, {
   3.8055025944080496`*^9, 3.805502624323163*^9}, {3.8055026772153587`*^9, 
   3.805502783940276*^9}, {3.805502854239046*^9, 3.80550289141193*^9}, 
   3.805539092142848*^9, 3.8055392008147535`*^9, {3.805539266359166*^9, 
   3.8055392950748816`*^9}, 3.8055476424770594`*^9, {3.805547743548753*^9, 
   3.805547782692649*^9}, 3.806061311121106*^9, {3.8060613762090282`*^9, 
   3.806061399563547*^9}, 3.806061485787599*^9, {3.8060618969884276`*^9, 
   3.8060619413353577`*^9}, 3.8060619883466463`*^9, 3.806062076635193*^9, {
   3.806062121638563*^9, 3.8060622477261877`*^9}, 3.8061436402582865`*^9, 
   3.812479908939375*^9, {3.8124805431737375`*^9, 3.8124805521521053`*^9}, 
   3.8127080260794463`*^9, 3.8127082080786448`*^9, 3.8127082767170734`*^9, 
   3.8127085542792435`*^9, 3.8127103778091154`*^9, 3.81271041590324*^9, 
   3.8127105663518934`*^9, 3.8127106029664683`*^9, 3.8127110945911274`*^9, {
   3.8127116204864073`*^9, 3.8127116438698387`*^9}, 3.8127120512537746`*^9, 
   3.8127132694632845`*^9, 3.812713528220541*^9},
 CellLabel->"Out[59]=",ExpressionUUID->"4c3d9b50-8f79-4e1e-9069-6b9d580a9863"],

Cell[BoxData["\<\"Expected Gain in SNB from regulating under the \
standard\"\>"], "Output",
 CellChangeTimes->{
  3.8054797335618596`*^9, 3.805479768529723*^9, 3.805479962547964*^9, 
   3.805480011263193*^9, 3.805480158628708*^9, 3.8054805207302823`*^9, {
   3.8054806478378954`*^9, 3.8054807227787147`*^9}, {3.805480754776679*^9, 
   3.805480802785015*^9}, {3.805480860379613*^9, 3.805480975306303*^9}, {
   3.805481055037532*^9, 3.805481138991852*^9}, {3.805481207897601*^9, 
   3.8054813233177137`*^9}, {3.805481668996467*^9, 3.8054816937761745`*^9}, {
   3.8054824596940994`*^9, 3.805482511675572*^9}, 3.80548270058698*^9, {
   3.805482938776911*^9, 3.805482957449237*^9}, 3.8054830084992285`*^9, 
   3.8055023938153105`*^9, 3.805502465059903*^9, 3.8055025620124936`*^9, {
   3.8055025944080496`*^9, 3.805502624323163*^9}, {3.8055026772153587`*^9, 
   3.805502783940276*^9}, {3.805502854239046*^9, 3.80550289141193*^9}, 
   3.805539092142848*^9, 3.8055392008147535`*^9, {3.805539266359166*^9, 
   3.8055392950748816`*^9}, 3.8055476424770594`*^9, {3.805547743548753*^9, 
   3.805547782692649*^9}, 3.806061311121106*^9, {3.8060613762090282`*^9, 
   3.806061399563547*^9}, 3.806061485787599*^9, {3.8060618969884276`*^9, 
   3.8060619413353577`*^9}, 3.8060619883466463`*^9, 3.806062076635193*^9, {
   3.806062121638563*^9, 3.8060622477261877`*^9}, 3.8061436402582865`*^9, 
   3.812479908939375*^9, {3.8124805431737375`*^9, 3.8124805521521053`*^9}, 
   3.8127080260794463`*^9, 3.8127082080786448`*^9, 3.8127082767170734`*^9, 
   3.8127085542792435`*^9, 3.8127103778091154`*^9, 3.81271041590324*^9, 
   3.8127105663518934`*^9, 3.8127106029664683`*^9, 3.8127110945911274`*^9, {
   3.8127116204864073`*^9, 3.8127116438698387`*^9}, 3.8127120512537746`*^9, 
   3.8127132694632845`*^9, 3.8127135282235327`*^9},
 CellLabel->"Out[60]=",ExpressionUUID->"340e7b4c-8838-436f-861c-9d191d0a75ae"],

Cell[BoxData["18188.476353872862`"], "Output",
 CellChangeTimes->{
  3.8054797335618596`*^9, 3.805479768529723*^9, 3.805479962547964*^9, 
   3.805480011263193*^9, 3.805480158628708*^9, 3.8054805207302823`*^9, {
   3.8054806478378954`*^9, 3.8054807227787147`*^9}, {3.805480754776679*^9, 
   3.805480802785015*^9}, {3.805480860379613*^9, 3.805480975306303*^9}, {
   3.805481055037532*^9, 3.805481138991852*^9}, {3.805481207897601*^9, 
   3.8054813233177137`*^9}, {3.805481668996467*^9, 3.8054816937761745`*^9}, {
   3.8054824596940994`*^9, 3.805482511675572*^9}, 3.80548270058698*^9, {
   3.805482938776911*^9, 3.805482957449237*^9}, 3.8054830084992285`*^9, 
   3.8055023938153105`*^9, 3.805502465059903*^9, 3.8055025620124936`*^9, {
   3.8055025944080496`*^9, 3.805502624323163*^9}, {3.8055026772153587`*^9, 
   3.805502783940276*^9}, {3.805502854239046*^9, 3.80550289141193*^9}, 
   3.805539092142848*^9, 3.8055392008147535`*^9, {3.805539266359166*^9, 
   3.8055392950748816`*^9}, 3.8055476424770594`*^9, {3.805547743548753*^9, 
   3.805547782692649*^9}, 3.806061311121106*^9, {3.8060613762090282`*^9, 
   3.806061399563547*^9}, 3.806061485787599*^9, {3.8060618969884276`*^9, 
   3.8060619413353577`*^9}, 3.8060619883466463`*^9, 3.806062076635193*^9, {
   3.806062121638563*^9, 3.8060622477261877`*^9}, 3.8061436402582865`*^9, 
   3.812479908939375*^9, {3.8124805431737375`*^9, 3.8124805521521053`*^9}, 
   3.8127080260794463`*^9, 3.8127082080786448`*^9, 3.8127082767170734`*^9, 
   3.8127085542792435`*^9, 3.8127103778091154`*^9, 3.81271041590324*^9, 
   3.8127105663518934`*^9, 3.8127106029664683`*^9, 3.8127110945911274`*^9, {
   3.8127116204864073`*^9, 3.8127116438698387`*^9}, 3.8127120512537746`*^9, 
   3.8127132694632845`*^9, 3.8127135282255273`*^9},
 CellLabel->"Out[61]=",ExpressionUUID->"c1a82b20-1a94-423e-85bb-9457ee3bb4d2"],

Cell[BoxData["\<\"Covariance between \[Alpha] and \[Gamma]\"\>"], "Output",
 CellChangeTimes->{
  3.8054797335618596`*^9, 3.805479768529723*^9, 3.805479962547964*^9, 
   3.805480011263193*^9, 3.805480158628708*^9, 3.8054805207302823`*^9, {
   3.8054806478378954`*^9, 3.8054807227787147`*^9}, {3.805480754776679*^9, 
   3.805480802785015*^9}, {3.805480860379613*^9, 3.805480975306303*^9}, {
   3.805481055037532*^9, 3.805481138991852*^9}, {3.805481207897601*^9, 
   3.8054813233177137`*^9}, {3.805481668996467*^9, 3.8054816937761745`*^9}, {
   3.8054824596940994`*^9, 3.805482511675572*^9}, 3.80548270058698*^9, {
   3.805482938776911*^9, 3.805482957449237*^9}, 3.8054830084992285`*^9, 
   3.8055023938153105`*^9, 3.805502465059903*^9, 3.8055025620124936`*^9, {
   3.8055025944080496`*^9, 3.805502624323163*^9}, {3.8055026772153587`*^9, 
   3.805502783940276*^9}, {3.805502854239046*^9, 3.80550289141193*^9}, 
   3.805539092142848*^9, 3.8055392008147535`*^9, {3.805539266359166*^9, 
   3.8055392950748816`*^9}, 3.8055476424770594`*^9, {3.805547743548753*^9, 
   3.805547782692649*^9}, 3.806061311121106*^9, {3.8060613762090282`*^9, 
   3.806061399563547*^9}, 3.806061485787599*^9, {3.8060618969884276`*^9, 
   3.8060619413353577`*^9}, 3.8060619883466463`*^9, 3.806062076635193*^9, {
   3.806062121638563*^9, 3.8060622477261877`*^9}, 3.8061436402582865`*^9, 
   3.812479908939375*^9, {3.8124805431737375`*^9, 3.8124805521521053`*^9}, 
   3.8127080260794463`*^9, 3.8127082080786448`*^9, 3.8127082767170734`*^9, 
   3.8127085542792435`*^9, 3.8127103778091154`*^9, 3.81271041590324*^9, 
   3.8127105663518934`*^9, 3.8127106029664683`*^9, 3.8127110945911274`*^9, {
   3.8127116204864073`*^9, 3.8127116438698387`*^9}, 3.8127120512537746`*^9, 
   3.8127132694632845`*^9, 3.812713528227522*^9},
 CellLabel->"Out[62]=",ExpressionUUID->"bf40d956-5cfa-42ba-bb08-50a28e3c4a83"],

Cell[BoxData[
 RowBox[{"-", "142.10526315789468`"}]], "Output",
 CellChangeTimes->{
  3.8054797335618596`*^9, 3.805479768529723*^9, 3.805479962547964*^9, 
   3.805480011263193*^9, 3.805480158628708*^9, 3.8054805207302823`*^9, {
   3.8054806478378954`*^9, 3.8054807227787147`*^9}, {3.805480754776679*^9, 
   3.805480802785015*^9}, {3.805480860379613*^9, 3.805480975306303*^9}, {
   3.805481055037532*^9, 3.805481138991852*^9}, {3.805481207897601*^9, 
   3.8054813233177137`*^9}, {3.805481668996467*^9, 3.8054816937761745`*^9}, {
   3.8054824596940994`*^9, 3.805482511675572*^9}, 3.80548270058698*^9, {
   3.805482938776911*^9, 3.805482957449237*^9}, 3.8054830084992285`*^9, 
   3.8055023938153105`*^9, 3.805502465059903*^9, 3.8055025620124936`*^9, {
   3.8055025944080496`*^9, 3.805502624323163*^9}, {3.8055026772153587`*^9, 
   3.805502783940276*^9}, {3.805502854239046*^9, 3.80550289141193*^9}, 
   3.805539092142848*^9, 3.8055392008147535`*^9, {3.805539266359166*^9, 
   3.8055392950748816`*^9}, 3.8055476424770594`*^9, {3.805547743548753*^9, 
   3.805547782692649*^9}, 3.806061311121106*^9, {3.8060613762090282`*^9, 
   3.806061399563547*^9}, 3.806061485787599*^9, {3.8060618969884276`*^9, 
   3.8060619413353577`*^9}, 3.8060619883466463`*^9, 3.806062076635193*^9, {
   3.806062121638563*^9, 3.8060622477261877`*^9}, 3.8061436402582865`*^9, 
   3.812479908939375*^9, {3.8124805431737375`*^9, 3.8124805521521053`*^9}, 
   3.8127080260794463`*^9, 3.8127082080786448`*^9, 3.8127082767170734`*^9, 
   3.8127085542792435`*^9, 3.8127103778091154`*^9, 3.81271041590324*^9, 
   3.8127105663518934`*^9, 3.8127106029664683`*^9, 3.8127110945911274`*^9, {
   3.8127116204864073`*^9, 3.8127116438698387`*^9}, 3.8127120512537746`*^9, 
   3.8127132694632845`*^9, 3.812713528230514*^9},
 CellLabel->"Out[63]=",ExpressionUUID->"9708c0bb-d6ce-4ee2-bac0-0cbf909b0013"],

Cell[BoxData["\<\"true tax\"\>"], "Output",
 CellChangeTimes->{
  3.8054797335618596`*^9, 3.805479768529723*^9, 3.805479962547964*^9, 
   3.805480011263193*^9, 3.805480158628708*^9, 3.8054805207302823`*^9, {
   3.8054806478378954`*^9, 3.8054807227787147`*^9}, {3.805480754776679*^9, 
   3.805480802785015*^9}, {3.805480860379613*^9, 3.805480975306303*^9}, {
   3.805481055037532*^9, 3.805481138991852*^9}, {3.805481207897601*^9, 
   3.8054813233177137`*^9}, {3.805481668996467*^9, 3.8054816937761745`*^9}, {
   3.8054824596940994`*^9, 3.805482511675572*^9}, 3.80548270058698*^9, {
   3.805482938776911*^9, 3.805482957449237*^9}, 3.8054830084992285`*^9, 
   3.8055023938153105`*^9, 3.805502465059903*^9, 3.8055025620124936`*^9, {
   3.8055025944080496`*^9, 3.805502624323163*^9}, {3.8055026772153587`*^9, 
   3.805502783940276*^9}, {3.805502854239046*^9, 3.80550289141193*^9}, 
   3.805539092142848*^9, 3.8055392008147535`*^9, {3.805539266359166*^9, 
   3.8055392950748816`*^9}, 3.8055476424770594`*^9, {3.805547743548753*^9, 
   3.805547782692649*^9}, 3.806061311121106*^9, {3.8060613762090282`*^9, 
   3.806061399563547*^9}, 3.806061485787599*^9, {3.8060618969884276`*^9, 
   3.8060619413353577`*^9}, 3.8060619883466463`*^9, 3.806062076635193*^9, {
   3.806062121638563*^9, 3.8060622477261877`*^9}, 3.8061436402582865`*^9, 
   3.812479908939375*^9, {3.8124805431737375`*^9, 3.8124805521521053`*^9}, 
   3.8127080260794463`*^9, 3.8127082080786448`*^9, 3.8127082767170734`*^9, 
   3.8127085542792435`*^9, 3.8127103778091154`*^9, 3.81271041590324*^9, 
   3.8127105663518934`*^9, 3.8127106029664683`*^9, 3.8127110945911274`*^9, {
   3.8127116204864073`*^9, 3.8127116438698387`*^9}, 3.8127120512537746`*^9, 
   3.8127132694632845`*^9, 3.8127135282325087`*^9},
 CellLabel->"Out[64]=",ExpressionUUID->"720e5528-627c-46be-bdb6-ff14cec00545"],

Cell[BoxData["743.3090464791178`"], "Output",
 CellChangeTimes->{
  3.8054797335618596`*^9, 3.805479768529723*^9, 3.805479962547964*^9, 
   3.805480011263193*^9, 3.805480158628708*^9, 3.8054805207302823`*^9, {
   3.8054806478378954`*^9, 3.8054807227787147`*^9}, {3.805480754776679*^9, 
   3.805480802785015*^9}, {3.805480860379613*^9, 3.805480975306303*^9}, {
   3.805481055037532*^9, 3.805481138991852*^9}, {3.805481207897601*^9, 
   3.8054813233177137`*^9}, {3.805481668996467*^9, 3.8054816937761745`*^9}, {
   3.8054824596940994`*^9, 3.805482511675572*^9}, 3.80548270058698*^9, {
   3.805482938776911*^9, 3.805482957449237*^9}, 3.8054830084992285`*^9, 
   3.8055023938153105`*^9, 3.805502465059903*^9, 3.8055025620124936`*^9, {
   3.8055025944080496`*^9, 3.805502624323163*^9}, {3.8055026772153587`*^9, 
   3.805502783940276*^9}, {3.805502854239046*^9, 3.80550289141193*^9}, 
   3.805539092142848*^9, 3.8055392008147535`*^9, {3.805539266359166*^9, 
   3.8055392950748816`*^9}, 3.8055476424770594`*^9, {3.805547743548753*^9, 
   3.805547782692649*^9}, 3.806061311121106*^9, {3.8060613762090282`*^9, 
   3.806061399563547*^9}, 3.806061485787599*^9, {3.8060618969884276`*^9, 
   3.8060619413353577`*^9}, 3.8060619883466463`*^9, 3.806062076635193*^9, {
   3.806062121638563*^9, 3.8060622477261877`*^9}, 3.8061436402582865`*^9, 
   3.812479908939375*^9, {3.8124805431737375`*^9, 3.8124805521521053`*^9}, 
   3.8127080260794463`*^9, 3.8127082080786448`*^9, 3.8127082767170734`*^9, 
   3.8127085542792435`*^9, 3.8127103778091154`*^9, 3.81271041590324*^9, 
   3.8127105663518934`*^9, 3.8127106029664683`*^9, 3.8127110945911274`*^9, {
   3.8127116204864073`*^9, 3.8127116438698387`*^9}, 3.8127120512537746`*^9, 
   3.8127132694632845`*^9, 3.8127135282345033`*^9},
 CellLabel->"Out[65]=",ExpressionUUID->"1472e81b-1498-4d04-bfbb-a33feedcf84a"],

Cell[BoxData["\<\"r\"\>"], "Output",
 CellChangeTimes->{
  3.8054797335618596`*^9, 3.805479768529723*^9, 3.805479962547964*^9, 
   3.805480011263193*^9, 3.805480158628708*^9, 3.8054805207302823`*^9, {
   3.8054806478378954`*^9, 3.8054807227787147`*^9}, {3.805480754776679*^9, 
   3.805480802785015*^9}, {3.805480860379613*^9, 3.805480975306303*^9}, {
   3.805481055037532*^9, 3.805481138991852*^9}, {3.805481207897601*^9, 
   3.8054813233177137`*^9}, {3.805481668996467*^9, 3.8054816937761745`*^9}, {
   3.8054824596940994`*^9, 3.805482511675572*^9}, 3.80548270058698*^9, {
   3.805482938776911*^9, 3.805482957449237*^9}, 3.8054830084992285`*^9, 
   3.8055023938153105`*^9, 3.805502465059903*^9, 3.8055025620124936`*^9, {
   3.8055025944080496`*^9, 3.805502624323163*^9}, {3.8055026772153587`*^9, 
   3.805502783940276*^9}, {3.805502854239046*^9, 3.80550289141193*^9}, 
   3.805539092142848*^9, 3.8055392008147535`*^9, {3.805539266359166*^9, 
   3.8055392950748816`*^9}, 3.8055476424770594`*^9, {3.805547743548753*^9, 
   3.805547782692649*^9}, 3.806061311121106*^9, {3.8060613762090282`*^9, 
   3.806061399563547*^9}, 3.806061485787599*^9, {3.8060618969884276`*^9, 
   3.8060619413353577`*^9}, 3.8060619883466463`*^9, 3.806062076635193*^9, {
   3.806062121638563*^9, 3.8060622477261877`*^9}, 3.8061436402582865`*^9, 
   3.812479908939375*^9, {3.8124805431737375`*^9, 3.8124805521521053`*^9}, 
   3.8127080260794463`*^9, 3.8127082080786448`*^9, 3.8127082767170734`*^9, 
   3.8127085542792435`*^9, 3.8127103778091154`*^9, 3.81271041590324*^9, 
   3.8127105663518934`*^9, 3.8127106029664683`*^9, 3.8127110945911274`*^9, {
   3.8127116204864073`*^9, 3.8127116438698387`*^9}, 3.8127120512537746`*^9, 
   3.8127132694632845`*^9, 3.812713528236498*^9},
 CellLabel->"Out[66]=",ExpressionUUID->"c93e01df-33e2-40c3-8193-a3bc52f2831b"],

Cell[BoxData["269.62483595280196`"], "Output",
 CellChangeTimes->{
  3.8054797335618596`*^9, 3.805479768529723*^9, 3.805479962547964*^9, 
   3.805480011263193*^9, 3.805480158628708*^9, 3.8054805207302823`*^9, {
   3.8054806478378954`*^9, 3.8054807227787147`*^9}, {3.805480754776679*^9, 
   3.805480802785015*^9}, {3.805480860379613*^9, 3.805480975306303*^9}, {
   3.805481055037532*^9, 3.805481138991852*^9}, {3.805481207897601*^9, 
   3.8054813233177137`*^9}, {3.805481668996467*^9, 3.8054816937761745`*^9}, {
   3.8054824596940994`*^9, 3.805482511675572*^9}, 3.80548270058698*^9, {
   3.805482938776911*^9, 3.805482957449237*^9}, 3.8054830084992285`*^9, 
   3.8055023938153105`*^9, 3.805502465059903*^9, 3.8055025620124936`*^9, {
   3.8055025944080496`*^9, 3.805502624323163*^9}, {3.8055026772153587`*^9, 
   3.805502783940276*^9}, {3.805502854239046*^9, 3.80550289141193*^9}, 
   3.805539092142848*^9, 3.8055392008147535`*^9, {3.805539266359166*^9, 
   3.8055392950748816`*^9}, 3.8055476424770594`*^9, {3.805547743548753*^9, 
   3.805547782692649*^9}, 3.806061311121106*^9, {3.8060613762090282`*^9, 
   3.806061399563547*^9}, 3.806061485787599*^9, {3.8060618969884276`*^9, 
   3.8060619413353577`*^9}, 3.8060619883466463`*^9, 3.806062076635193*^9, {
   3.806062121638563*^9, 3.8060622477261877`*^9}, 3.8061436402582865`*^9, 
   3.812479908939375*^9, {3.8124805431737375`*^9, 3.8124805521521053`*^9}, 
   3.8127080260794463`*^9, 3.8127082080786448`*^9, 3.8127082767170734`*^9, 
   3.8127085542792435`*^9, 3.8127103778091154`*^9, 3.81271041590324*^9, 
   3.8127105663518934`*^9, 3.8127106029664683`*^9, 3.8127110945911274`*^9, {
   3.8127116204864073`*^9, 3.8127116438698387`*^9}, 3.8127120512537746`*^9, 
   3.8127132694632845`*^9, 3.8127135282394896`*^9},
 CellLabel->"Out[67]=",ExpressionUUID->"1877f085-bd5b-4a44-869f-43c27c956ffe"],

Cell[BoxData["\<\"ESNB True \[Tau]\"\>"], "Output",
 CellChangeTimes->{
  3.8054797335618596`*^9, 3.805479768529723*^9, 3.805479962547964*^9, 
   3.805480011263193*^9, 3.805480158628708*^9, 3.8054805207302823`*^9, {
   3.8054806478378954`*^9, 3.8054807227787147`*^9}, {3.805480754776679*^9, 
   3.805480802785015*^9}, {3.805480860379613*^9, 3.805480975306303*^9}, {
   3.805481055037532*^9, 3.805481138991852*^9}, {3.805481207897601*^9, 
   3.8054813233177137`*^9}, {3.805481668996467*^9, 3.8054816937761745`*^9}, {
   3.8054824596940994`*^9, 3.805482511675572*^9}, 3.80548270058698*^9, {
   3.805482938776911*^9, 3.805482957449237*^9}, 3.8054830084992285`*^9, 
   3.8055023938153105`*^9, 3.805502465059903*^9, 3.8055025620124936`*^9, {
   3.8055025944080496`*^9, 3.805502624323163*^9}, {3.8055026772153587`*^9, 
   3.805502783940276*^9}, {3.805502854239046*^9, 3.80550289141193*^9}, 
   3.805539092142848*^9, 3.8055392008147535`*^9, {3.805539266359166*^9, 
   3.8055392950748816`*^9}, 3.8055476424770594`*^9, {3.805547743548753*^9, 
   3.805547782692649*^9}, 3.806061311121106*^9, {3.8060613762090282`*^9, 
   3.806061399563547*^9}, 3.806061485787599*^9, {3.8060618969884276`*^9, 
   3.8060619413353577`*^9}, 3.8060619883466463`*^9, 3.806062076635193*^9, {
   3.806062121638563*^9, 3.8060622477261877`*^9}, 3.8061436402582865`*^9, 
   3.812479908939375*^9, {3.8124805431737375`*^9, 3.8124805521521053`*^9}, 
   3.8127080260794463`*^9, 3.8127082080786448`*^9, 3.8127082767170734`*^9, 
   3.8127085542792435`*^9, 3.8127103778091154`*^9, 3.81271041590324*^9, 
   3.8127105663518934`*^9, 3.8127106029664683`*^9, 3.8127110945911274`*^9, {
   3.8127116204864073`*^9, 3.8127116438698387`*^9}, 3.8127120512537746`*^9, 
   3.8127132694632845`*^9, 3.8127135282414846`*^9},
 CellLabel->"Out[68]=",ExpressionUUID->"968c2c9c-142f-4115-8c3f-423088b70b44"],

Cell[BoxData["2047.7905364220262`"], "Output",
 CellChangeTimes->{
  3.8054797335618596`*^9, 3.805479768529723*^9, 3.805479962547964*^9, 
   3.805480011263193*^9, 3.805480158628708*^9, 3.8054805207302823`*^9, {
   3.8054806478378954`*^9, 3.8054807227787147`*^9}, {3.805480754776679*^9, 
   3.805480802785015*^9}, {3.805480860379613*^9, 3.805480975306303*^9}, {
   3.805481055037532*^9, 3.805481138991852*^9}, {3.805481207897601*^9, 
   3.8054813233177137`*^9}, {3.805481668996467*^9, 3.8054816937761745`*^9}, {
   3.8054824596940994`*^9, 3.805482511675572*^9}, 3.80548270058698*^9, {
   3.805482938776911*^9, 3.805482957449237*^9}, 3.8054830084992285`*^9, 
   3.8055023938153105`*^9, 3.805502465059903*^9, 3.8055025620124936`*^9, {
   3.8055025944080496`*^9, 3.805502624323163*^9}, {3.8055026772153587`*^9, 
   3.805502783940276*^9}, {3.805502854239046*^9, 3.80550289141193*^9}, 
   3.805539092142848*^9, 3.8055392008147535`*^9, {3.805539266359166*^9, 
   3.8055392950748816`*^9}, 3.8055476424770594`*^9, {3.805547743548753*^9, 
   3.805547782692649*^9}, 3.806061311121106*^9, {3.8060613762090282`*^9, 
   3.806061399563547*^9}, 3.806061485787599*^9, {3.8060618969884276`*^9, 
   3.8060619413353577`*^9}, 3.8060619883466463`*^9, 3.806062076635193*^9, {
   3.806062121638563*^9, 3.8060622477261877`*^9}, 3.8061436402582865`*^9, 
   3.812479908939375*^9, {3.8124805431737375`*^9, 3.8124805521521053`*^9}, 
   3.8127080260794463`*^9, 3.8127082080786448`*^9, 3.8127082767170734`*^9, 
   3.8127085542792435`*^9, 3.8127103778091154`*^9, 3.81271041590324*^9, 
   3.8127105663518934`*^9, 3.8127106029664683`*^9, 3.8127110945911274`*^9, {
   3.8127116204864073`*^9, 3.8127116438698387`*^9}, 3.8127120512537746`*^9, 
   3.8127132694632845`*^9, 3.8127135282434793`*^9},
 CellLabel->"Out[69]=",ExpressionUUID->"b16690f5-ed57-4a36-9dd8-e24347c95e7f"],

Cell[BoxData["\<\"Expected Gain in SNB from regulating under the true \
tax\"\>"], "Output",
 CellChangeTimes->{
  3.8054797335618596`*^9, 3.805479768529723*^9, 3.805479962547964*^9, 
   3.805480011263193*^9, 3.805480158628708*^9, 3.8054805207302823`*^9, {
   3.8054806478378954`*^9, 3.8054807227787147`*^9}, {3.805480754776679*^9, 
   3.805480802785015*^9}, {3.805480860379613*^9, 3.805480975306303*^9}, {
   3.805481055037532*^9, 3.805481138991852*^9}, {3.805481207897601*^9, 
   3.8054813233177137`*^9}, {3.805481668996467*^9, 3.8054816937761745`*^9}, {
   3.8054824596940994`*^9, 3.805482511675572*^9}, 3.80548270058698*^9, {
   3.805482938776911*^9, 3.805482957449237*^9}, 3.8054830084992285`*^9, 
   3.8055023938153105`*^9, 3.805502465059903*^9, 3.8055025620124936`*^9, {
   3.8055025944080496`*^9, 3.805502624323163*^9}, {3.8055026772153587`*^9, 
   3.805502783940276*^9}, {3.805502854239046*^9, 3.80550289141193*^9}, 
   3.805539092142848*^9, 3.8055392008147535`*^9, {3.805539266359166*^9, 
   3.8055392950748816`*^9}, 3.8055476424770594`*^9, {3.805547743548753*^9, 
   3.805547782692649*^9}, 3.806061311121106*^9, {3.8060613762090282`*^9, 
   3.806061399563547*^9}, 3.806061485787599*^9, {3.8060618969884276`*^9, 
   3.8060619413353577`*^9}, 3.8060619883466463`*^9, 3.806062076635193*^9, {
   3.806062121638563*^9, 3.8060622477261877`*^9}, 3.8061436402582865`*^9, 
   3.812479908939375*^9, {3.8124805431737375`*^9, 3.8124805521521053`*^9}, 
   3.8127080260794463`*^9, 3.8127082080786448`*^9, 3.8127082767170734`*^9, 
   3.8127085542792435`*^9, 3.8127103778091154`*^9, 3.81271041590324*^9, 
   3.8127105663518934`*^9, 3.8127106029664683`*^9, 3.8127110945911274`*^9, {
   3.8127116204864073`*^9, 3.8127116438698387`*^9}, 3.8127120512537746`*^9, 
   3.8127132694632845`*^9, 3.8127135282464714`*^9},
 CellLabel->"Out[70]=",ExpressionUUID->"e6cd8845-c0a8-48b3-a185-cb698677d5b8"],

Cell[BoxData["17604.68794292647`"], "Output",
 CellChangeTimes->{
  3.8054797335618596`*^9, 3.805479768529723*^9, 3.805479962547964*^9, 
   3.805480011263193*^9, 3.805480158628708*^9, 3.8054805207302823`*^9, {
   3.8054806478378954`*^9, 3.8054807227787147`*^9}, {3.805480754776679*^9, 
   3.805480802785015*^9}, {3.805480860379613*^9, 3.805480975306303*^9}, {
   3.805481055037532*^9, 3.805481138991852*^9}, {3.805481207897601*^9, 
   3.8054813233177137`*^9}, {3.805481668996467*^9, 3.8054816937761745`*^9}, {
   3.8054824596940994`*^9, 3.805482511675572*^9}, 3.80548270058698*^9, {
   3.805482938776911*^9, 3.805482957449237*^9}, 3.8054830084992285`*^9, 
   3.8055023938153105`*^9, 3.805502465059903*^9, 3.8055025620124936`*^9, {
   3.8055025944080496`*^9, 3.805502624323163*^9}, {3.8055026772153587`*^9, 
   3.805502783940276*^9}, {3.805502854239046*^9, 3.80550289141193*^9}, 
   3.805539092142848*^9, 3.8055392008147535`*^9, {3.805539266359166*^9, 
   3.8055392950748816`*^9}, 3.8055476424770594`*^9, {3.805547743548753*^9, 
   3.805547782692649*^9}, 3.806061311121106*^9, {3.8060613762090282`*^9, 
   3.806061399563547*^9}, 3.806061485787599*^9, {3.8060618969884276`*^9, 
   3.8060619413353577`*^9}, 3.8060619883466463`*^9, 3.806062076635193*^9, {
   3.806062121638563*^9, 3.8060622477261877`*^9}, 3.8061436402582865`*^9, 
   3.812479908939375*^9, {3.8124805431737375`*^9, 3.8124805521521053`*^9}, 
   3.8127080260794463`*^9, 3.8127082080786448`*^9, 3.8127082767170734`*^9, 
   3.8127085542792435`*^9, 3.8127103778091154`*^9, 3.81271041590324*^9, 
   3.8127105663518934`*^9, 3.8127106029664683`*^9, 3.8127110945911274`*^9, {
   3.8127116204864073`*^9, 3.8127116438698387`*^9}, 3.8127120512537746`*^9, 
   3.8127132694632845`*^9, 3.8127135282484655`*^9},
 CellLabel->"Out[71]=",ExpressionUUID->"6f901064-f319-4180-939b-28c74dee252f"],

Cell[BoxData["\<\"True \[CapitalDelta]\"\>"], "Output",
 CellChangeTimes->{
  3.8054797335618596`*^9, 3.805479768529723*^9, 3.805479962547964*^9, 
   3.805480011263193*^9, 3.805480158628708*^9, 3.8054805207302823`*^9, {
   3.8054806478378954`*^9, 3.8054807227787147`*^9}, {3.805480754776679*^9, 
   3.805480802785015*^9}, {3.805480860379613*^9, 3.805480975306303*^9}, {
   3.805481055037532*^9, 3.805481138991852*^9}, {3.805481207897601*^9, 
   3.8054813233177137`*^9}, {3.805481668996467*^9, 3.8054816937761745`*^9}, {
   3.8054824596940994`*^9, 3.805482511675572*^9}, 3.80548270058698*^9, {
   3.805482938776911*^9, 3.805482957449237*^9}, 3.8054830084992285`*^9, 
   3.8055023938153105`*^9, 3.805502465059903*^9, 3.8055025620124936`*^9, {
   3.8055025944080496`*^9, 3.805502624323163*^9}, {3.8055026772153587`*^9, 
   3.805502783940276*^9}, {3.805502854239046*^9, 3.80550289141193*^9}, 
   3.805539092142848*^9, 3.8055392008147535`*^9, {3.805539266359166*^9, 
   3.8055392950748816`*^9}, 3.8055476424770594`*^9, {3.805547743548753*^9, 
   3.805547782692649*^9}, 3.806061311121106*^9, {3.8060613762090282`*^9, 
   3.806061399563547*^9}, 3.806061485787599*^9, {3.8060618969884276`*^9, 
   3.8060619413353577`*^9}, 3.8060619883466463`*^9, 3.806062076635193*^9, {
   3.806062121638563*^9, 3.8060622477261877`*^9}, 3.8061436402582865`*^9, 
   3.812479908939375*^9, {3.8124805431737375`*^9, 3.8124805521521053`*^9}, 
   3.8127080260794463`*^9, 3.8127082080786448`*^9, 3.8127082767170734`*^9, 
   3.8127085542792435`*^9, 3.8127103778091154`*^9, 3.81271041590324*^9, 
   3.8127105663518934`*^9, 3.8127106029664683`*^9, 3.8127110945911274`*^9, {
   3.8127116204864073`*^9, 3.8127116438698387`*^9}, 3.8127120512537746`*^9, 
   3.8127132694632845`*^9, 3.8127135282504606`*^9},
 CellLabel->"Out[72]=",ExpressionUUID->"db9d170e-6b90-44d3-a61c-a2db2bd71a46"],

Cell[BoxData[
 RowBox[{"-", "583.7884109463946`"}]], "Output",
 CellChangeTimes->{
  3.8054797335618596`*^9, 3.805479768529723*^9, 3.805479962547964*^9, 
   3.805480011263193*^9, 3.805480158628708*^9, 3.8054805207302823`*^9, {
   3.8054806478378954`*^9, 3.8054807227787147`*^9}, {3.805480754776679*^9, 
   3.805480802785015*^9}, {3.805480860379613*^9, 3.805480975306303*^9}, {
   3.805481055037532*^9, 3.805481138991852*^9}, {3.805481207897601*^9, 
   3.8054813233177137`*^9}, {3.805481668996467*^9, 3.8054816937761745`*^9}, {
   3.8054824596940994`*^9, 3.805482511675572*^9}, 3.80548270058698*^9, {
   3.805482938776911*^9, 3.805482957449237*^9}, 3.8054830084992285`*^9, 
   3.8055023938153105`*^9, 3.805502465059903*^9, 3.8055025620124936`*^9, {
   3.8055025944080496`*^9, 3.805502624323163*^9}, {3.8055026772153587`*^9, 
   3.805502783940276*^9}, {3.805502854239046*^9, 3.80550289141193*^9}, 
   3.805539092142848*^9, 3.8055392008147535`*^9, {3.805539266359166*^9, 
   3.8055392950748816`*^9}, 3.8055476424770594`*^9, {3.805547743548753*^9, 
   3.805547782692649*^9}, 3.806061311121106*^9, {3.8060613762090282`*^9, 
   3.806061399563547*^9}, 3.806061485787599*^9, {3.8060618969884276`*^9, 
   3.8060619413353577`*^9}, 3.8060619883466463`*^9, 3.806062076635193*^9, {
   3.806062121638563*^9, 3.8060622477261877`*^9}, 3.8061436402582865`*^9, 
   3.812479908939375*^9, {3.8124805431737375`*^9, 3.8124805521521053`*^9}, 
   3.8127080260794463`*^9, 3.8127082080786448`*^9, 3.8127082767170734`*^9, 
   3.8127085542792435`*^9, 3.8127103778091154`*^9, 3.81271041590324*^9, 
   3.8127105663518934`*^9, 3.8127106029664683`*^9, 3.8127110945911274`*^9, {
   3.8127116204864073`*^9, 3.8127116438698387`*^9}, 3.8127120512537746`*^9, 
   3.8127132694632845`*^9, 3.8127135282534523`*^9},
 CellLabel->"Out[73]=",ExpressionUUID->"370ed5ee-b937-4bfa-8195-8d310424a36a"],

Cell[BoxData["\<\"Expected SNB for an unregulated system, assuming \[Gamma]=1\
\"\>"], "Output",
 CellChangeTimes->{
  3.8054797335618596`*^9, 3.805479768529723*^9, 3.805479962547964*^9, 
   3.805480011263193*^9, 3.805480158628708*^9, 3.8054805207302823`*^9, {
   3.8054806478378954`*^9, 3.8054807227787147`*^9}, {3.805480754776679*^9, 
   3.805480802785015*^9}, {3.805480860379613*^9, 3.805480975306303*^9}, {
   3.805481055037532*^9, 3.805481138991852*^9}, {3.805481207897601*^9, 
   3.8054813233177137`*^9}, {3.805481668996467*^9, 3.8054816937761745`*^9}, {
   3.8054824596940994`*^9, 3.805482511675572*^9}, 3.80548270058698*^9, {
   3.805482938776911*^9, 3.805482957449237*^9}, 3.8054830084992285`*^9, 
   3.8055023938153105`*^9, 3.805502465059903*^9, 3.8055025620124936`*^9, {
   3.8055025944080496`*^9, 3.805502624323163*^9}, {3.8055026772153587`*^9, 
   3.805502783940276*^9}, {3.805502854239046*^9, 3.80550289141193*^9}, 
   3.805539092142848*^9, 3.8055392008147535`*^9, {3.805539266359166*^9, 
   3.8055392950748816`*^9}, 3.8055476424770594`*^9, {3.805547743548753*^9, 
   3.805547782692649*^9}, 3.806061311121106*^9, {3.8060613762090282`*^9, 
   3.806061399563547*^9}, 3.806061485787599*^9, {3.8060618969884276`*^9, 
   3.8060619413353577`*^9}, 3.8060619883466463`*^9, 3.806062076635193*^9, {
   3.806062121638563*^9, 3.8060622477261877`*^9}, 3.8061436402582865`*^9, 
   3.812479908939375*^9, {3.8124805431737375`*^9, 3.8124805521521053`*^9}, 
   3.8127080260794463`*^9, 3.8127082080786448`*^9, 3.8127082767170734`*^9, 
   3.8127085542792435`*^9, 3.8127103778091154`*^9, 3.81271041590324*^9, 
   3.8127105663518934`*^9, 3.8127106029664683`*^9, 3.8127110945911274`*^9, {
   3.8127116204864073`*^9, 3.8127116438698387`*^9}, 3.8127120512537746`*^9, 
   3.8127132694632845`*^9, 3.812713528255447*^9},
 CellLabel->"Out[74]=",ExpressionUUID->"9b0d1a96-2b64-47a7-9db1-ab7c11bdda8c"],

Cell[BoxData["539.0627885503236`"], "Output",
 CellChangeTimes->{
  3.8054797335618596`*^9, 3.805479768529723*^9, 3.805479962547964*^9, 
   3.805480011263193*^9, 3.805480158628708*^9, 3.8054805207302823`*^9, {
   3.8054806478378954`*^9, 3.8054807227787147`*^9}, {3.805480754776679*^9, 
   3.805480802785015*^9}, {3.805480860379613*^9, 3.805480975306303*^9}, {
   3.805481055037532*^9, 3.805481138991852*^9}, {3.805481207897601*^9, 
   3.8054813233177137`*^9}, {3.805481668996467*^9, 3.8054816937761745`*^9}, {
   3.8054824596940994`*^9, 3.805482511675572*^9}, 3.80548270058698*^9, {
   3.805482938776911*^9, 3.805482957449237*^9}, 3.8054830084992285`*^9, 
   3.8055023938153105`*^9, 3.805502465059903*^9, 3.8055025620124936`*^9, {
   3.8055025944080496`*^9, 3.805502624323163*^9}, {3.8055026772153587`*^9, 
   3.805502783940276*^9}, {3.805502854239046*^9, 3.80550289141193*^9}, 
   3.805539092142848*^9, 3.8055392008147535`*^9, {3.805539266359166*^9, 
   3.8055392950748816`*^9}, 3.8055476424770594`*^9, {3.805547743548753*^9, 
   3.805547782692649*^9}, 3.806061311121106*^9, {3.8060613762090282`*^9, 
   3.806061399563547*^9}, 3.806061485787599*^9, {3.8060618969884276`*^9, 
   3.8060619413353577`*^9}, 3.8060619883466463`*^9, 3.806062076635193*^9, {
   3.806062121638563*^9, 3.8060622477261877`*^9}, 3.8061436402582865`*^9, 
   3.812479908939375*^9, {3.8124805431737375`*^9, 3.8124805521521053`*^9}, 
   3.8127080260794463`*^9, 3.8127082080786448`*^9, 3.8127082767170734`*^9, 
   3.8127085542792435`*^9, 3.8127103778091154`*^9, 3.81271041590324*^9, 
   3.8127105663518934`*^9, 3.8127106029664683`*^9, 3.8127110945911274`*^9, {
   3.8127116204864073`*^9, 3.8127116438698387`*^9}, 3.8127120512537746`*^9, 
   3.8127132694632845`*^9, 3.8127135282574415`*^9},
 CellLabel->"Out[75]=",ExpressionUUID->"ab5a2342-9bed-4ed7-9b85-c19365b0e2cf"],

Cell[BoxData["\<\"Weitzman's Tax = MD at ehat\"\>"], "Output",
 CellChangeTimes->{
  3.8054797335618596`*^9, 3.805479768529723*^9, 3.805479962547964*^9, 
   3.805480011263193*^9, 3.805480158628708*^9, 3.8054805207302823`*^9, {
   3.8054806478378954`*^9, 3.8054807227787147`*^9}, {3.805480754776679*^9, 
   3.805480802785015*^9}, {3.805480860379613*^9, 3.805480975306303*^9}, {
   3.805481055037532*^9, 3.805481138991852*^9}, {3.805481207897601*^9, 
   3.8054813233177137`*^9}, {3.805481668996467*^9, 3.8054816937761745`*^9}, {
   3.8054824596940994`*^9, 3.805482511675572*^9}, 3.80548270058698*^9, {
   3.805482938776911*^9, 3.805482957449237*^9}, 3.8054830084992285`*^9, 
   3.8055023938153105`*^9, 3.805502465059903*^9, 3.8055025620124936`*^9, {
   3.8055025944080496`*^9, 3.805502624323163*^9}, {3.8055026772153587`*^9, 
   3.805502783940276*^9}, {3.805502854239046*^9, 3.80550289141193*^9}, 
   3.805539092142848*^9, 3.8055392008147535`*^9, {3.805539266359166*^9, 
   3.8055392950748816`*^9}, 3.8055476424770594`*^9, {3.805547743548753*^9, 
   3.805547782692649*^9}, 3.806061311121106*^9, {3.8060613762090282`*^9, 
   3.806061399563547*^9}, 3.806061485787599*^9, {3.8060618969884276`*^9, 
   3.8060619413353577`*^9}, 3.8060619883466463`*^9, 3.806062076635193*^9, {
   3.806062121638563*^9, 3.8060622477261877`*^9}, 3.8061436402582865`*^9, 
   3.812479908939375*^9, {3.8124805431737375`*^9, 3.8124805521521053`*^9}, 
   3.8127080260794463`*^9, 3.8127082080786448`*^9, 3.8127082767170734`*^9, 
   3.8127085542792435`*^9, 3.8127103778091154`*^9, 3.81271041590324*^9, 
   3.8127105663518934`*^9, 3.8127106029664683`*^9, 3.8127110945911274`*^9, {
   3.8127116204864073`*^9, 3.8127116438698387`*^9}, 3.8127120512537746`*^9, 
   3.8127132694632845`*^9, 3.8127135282594366`*^9},
 CellLabel->"Out[76]=",ExpressionUUID->"9e13bdc1-29df-489e-8d09-ae1b3cf0d21b"],

Cell[BoxData["473.68421052631584`"], "Output",
 CellChangeTimes->{
  3.8054797335618596`*^9, 3.805479768529723*^9, 3.805479962547964*^9, 
   3.805480011263193*^9, 3.805480158628708*^9, 3.8054805207302823`*^9, {
   3.8054806478378954`*^9, 3.8054807227787147`*^9}, {3.805480754776679*^9, 
   3.805480802785015*^9}, {3.805480860379613*^9, 3.805480975306303*^9}, {
   3.805481055037532*^9, 3.805481138991852*^9}, {3.805481207897601*^9, 
   3.8054813233177137`*^9}, {3.805481668996467*^9, 3.8054816937761745`*^9}, {
   3.8054824596940994`*^9, 3.805482511675572*^9}, 3.80548270058698*^9, {
   3.805482938776911*^9, 3.805482957449237*^9}, 3.8054830084992285`*^9, 
   3.8055023938153105`*^9, 3.805502465059903*^9, 3.8055025620124936`*^9, {
   3.8055025944080496`*^9, 3.805502624323163*^9}, {3.8055026772153587`*^9, 
   3.805502783940276*^9}, {3.805502854239046*^9, 3.80550289141193*^9}, 
   3.805539092142848*^9, 3.8055392008147535`*^9, {3.805539266359166*^9, 
   3.8055392950748816`*^9}, 3.8055476424770594`*^9, {3.805547743548753*^9, 
   3.805547782692649*^9}, 3.806061311121106*^9, {3.8060613762090282`*^9, 
   3.806061399563547*^9}, 3.806061485787599*^9, {3.8060618969884276`*^9, 
   3.8060619413353577`*^9}, 3.8060619883466463`*^9, 3.806062076635193*^9, {
   3.806062121638563*^9, 3.8060622477261877`*^9}, 3.8061436402582865`*^9, 
   3.812479908939375*^9, {3.8124805431737375`*^9, 3.8124805521521053`*^9}, 
   3.8127080260794463`*^9, 3.8127082080786448`*^9, 3.8127082767170734`*^9, 
   3.8127085542792435`*^9, 3.8127103778091154`*^9, 3.81271041590324*^9, 
   3.8127105663518934`*^9, 3.8127106029664683`*^9, 3.8127110945911274`*^9, {
   3.8127116204864073`*^9, 3.8127116438698387`*^9}, 3.8127120512537746`*^9, 
   3.8127132694632845`*^9, 3.8127135282614307`*^9},
 CellLabel->"Out[77]=",ExpressionUUID->"122c57ef-ccb7-44e4-896a-f9cd3d9dd7aa"],

Cell[BoxData["\<\"ESNB Tax Weitman\"\>"], "Output",
 CellChangeTimes->{
  3.8054797335618596`*^9, 3.805479768529723*^9, 3.805479962547964*^9, 
   3.805480011263193*^9, 3.805480158628708*^9, 3.8054805207302823`*^9, {
   3.8054806478378954`*^9, 3.8054807227787147`*^9}, {3.805480754776679*^9, 
   3.805480802785015*^9}, {3.805480860379613*^9, 3.805480975306303*^9}, {
   3.805481055037532*^9, 3.805481138991852*^9}, {3.805481207897601*^9, 
   3.8054813233177137`*^9}, {3.805481668996467*^9, 3.8054816937761745`*^9}, {
   3.8054824596940994`*^9, 3.805482511675572*^9}, 3.80548270058698*^9, {
   3.805482938776911*^9, 3.805482957449237*^9}, 3.8054830084992285`*^9, 
   3.8055023938153105`*^9, 3.805502465059903*^9, 3.8055025620124936`*^9, {
   3.8055025944080496`*^9, 3.805502624323163*^9}, {3.8055026772153587`*^9, 
   3.805502783940276*^9}, {3.805502854239046*^9, 3.80550289141193*^9}, 
   3.805539092142848*^9, 3.8055392008147535`*^9, {3.805539266359166*^9, 
   3.8055392950748816`*^9}, 3.8055476424770594`*^9, {3.805547743548753*^9, 
   3.805547782692649*^9}, 3.806061311121106*^9, {3.8060613762090282`*^9, 
   3.806061399563547*^9}, 3.806061485787599*^9, {3.8060618969884276`*^9, 
   3.8060619413353577`*^9}, 3.8060619883466463`*^9, 3.806062076635193*^9, {
   3.806062121638563*^9, 3.8060622477261877`*^9}, 3.8061436402582865`*^9, 
   3.812479908939375*^9, {3.8124805431737375`*^9, 3.8124805521521053`*^9}, 
   3.8127080260794463`*^9, 3.8127082080786448`*^9, 3.8127082767170734`*^9, 
   3.8127085542792435`*^9, 3.8127103778091154`*^9, 3.81271041590324*^9, 
   3.8127105663518934`*^9, 3.8127106029664683`*^9, 3.8127110945911274`*^9, {
   3.8127116204864073`*^9, 3.8127116438698387`*^9}, 3.8127120512537746`*^9, 
   3.8127132694632845`*^9, 3.812713528264423*^9},
 CellLabel->"Out[78]=",ExpressionUUID->"db7abda0-f71c-41a2-9a62-ab41ce886a5e"],

Cell[BoxData["2670.6417359187444`"], "Output",
 CellChangeTimes->{
  3.8054797335618596`*^9, 3.805479768529723*^9, 3.805479962547964*^9, 
   3.805480011263193*^9, 3.805480158628708*^9, 3.8054805207302823`*^9, {
   3.8054806478378954`*^9, 3.8054807227787147`*^9}, {3.805480754776679*^9, 
   3.805480802785015*^9}, {3.805480860379613*^9, 3.805480975306303*^9}, {
   3.805481055037532*^9, 3.805481138991852*^9}, {3.805481207897601*^9, 
   3.8054813233177137`*^9}, {3.805481668996467*^9, 3.8054816937761745`*^9}, {
   3.8054824596940994`*^9, 3.805482511675572*^9}, 3.80548270058698*^9, {
   3.805482938776911*^9, 3.805482957449237*^9}, 3.8054830084992285`*^9, 
   3.8055023938153105`*^9, 3.805502465059903*^9, 3.8055025620124936`*^9, {
   3.8055025944080496`*^9, 3.805502624323163*^9}, {3.8055026772153587`*^9, 
   3.805502783940276*^9}, {3.805502854239046*^9, 3.80550289141193*^9}, 
   3.805539092142848*^9, 3.8055392008147535`*^9, {3.805539266359166*^9, 
   3.8055392950748816`*^9}, 3.8055476424770594`*^9, {3.805547743548753*^9, 
   3.805547782692649*^9}, 3.806061311121106*^9, {3.8060613762090282`*^9, 
   3.806061399563547*^9}, 3.806061485787599*^9, {3.8060618969884276`*^9, 
   3.8060619413353577`*^9}, 3.8060619883466463`*^9, 3.806062076635193*^9, {
   3.806062121638563*^9, 3.8060622477261877`*^9}, 3.8061436402582865`*^9, 
   3.812479908939375*^9, {3.8124805431737375`*^9, 3.8124805521521053`*^9}, 
   3.8127080260794463`*^9, 3.8127082080786448`*^9, 3.8127082767170734`*^9, 
   3.8127085542792435`*^9, 3.8127103778091154`*^9, 3.81271041590324*^9, 
   3.8127105663518934`*^9, 3.8127106029664683`*^9, 3.8127110945911274`*^9, {
   3.8127116204864073`*^9, 3.8127116438698387`*^9}, 3.8127120512537746`*^9, 
   3.8127132694632845`*^9, 3.812713528266431*^9},
 CellLabel->"Out[79]=",ExpressionUUID->"df710ae5-1710-4816-b52b-c1771bb5724f"],

Cell[BoxData["\<\"Expected Gain in SNB from regulating under Weitzman's tax \
and assuming \[Gamma]=1 when deriving unregulated ESNB\"\>"], "Output",
 CellChangeTimes->{
  3.8054797335618596`*^9, 3.805479768529723*^9, 3.805479962547964*^9, 
   3.805480011263193*^9, 3.805480158628708*^9, 3.8054805207302823`*^9, {
   3.8054806478378954`*^9, 3.8054807227787147`*^9}, {3.805480754776679*^9, 
   3.805480802785015*^9}, {3.805480860379613*^9, 3.805480975306303*^9}, {
   3.805481055037532*^9, 3.805481138991852*^9}, {3.805481207897601*^9, 
   3.8054813233177137`*^9}, {3.805481668996467*^9, 3.8054816937761745`*^9}, {
   3.8054824596940994`*^9, 3.805482511675572*^9}, 3.80548270058698*^9, {
   3.805482938776911*^9, 3.805482957449237*^9}, 3.8054830084992285`*^9, 
   3.8055023938153105`*^9, 3.805502465059903*^9, 3.8055025620124936`*^9, {
   3.8055025944080496`*^9, 3.805502624323163*^9}, {3.8055026772153587`*^9, 
   3.805502783940276*^9}, {3.805502854239046*^9, 3.80550289141193*^9}, 
   3.805539092142848*^9, 3.8055392008147535`*^9, {3.805539266359166*^9, 
   3.8055392950748816`*^9}, 3.8055476424770594`*^9, {3.805547743548753*^9, 
   3.805547782692649*^9}, 3.806061311121106*^9, {3.8060613762090282`*^9, 
   3.806061399563547*^9}, 3.806061485787599*^9, {3.8060618969884276`*^9, 
   3.8060619413353577`*^9}, 3.8060619883466463`*^9, 3.806062076635193*^9, {
   3.806062121638563*^9, 3.8060622477261877`*^9}, 3.8061436402582865`*^9, 
   3.812479908939375*^9, {3.8124805431737375`*^9, 3.8124805521521053`*^9}, 
   3.8127080260794463`*^9, 3.8127082080786448`*^9, 3.8127082767170734`*^9, 
   3.8127085542792435`*^9, 3.8127103778091154`*^9, 3.81271041590324*^9, 
   3.8127105663518934`*^9, 3.8127106029664683`*^9, 3.8127110945911274`*^9, {
   3.8127116204864073`*^9, 3.8127116438698387`*^9}, 3.8127120512537746`*^9, 
   3.8127132694632845`*^9, 3.812713528268412*^9},
 CellLabel->"Out[80]=",ExpressionUUID->"73d9a939-b8b8-4af2-be6b-ca643ccbdf2e"],

Cell[BoxData["2131.578947368421`"], "Output",
 CellChangeTimes->{
  3.8054797335618596`*^9, 3.805479768529723*^9, 3.805479962547964*^9, 
   3.805480011263193*^9, 3.805480158628708*^9, 3.8054805207302823`*^9, {
   3.8054806478378954`*^9, 3.8054807227787147`*^9}, {3.805480754776679*^9, 
   3.805480802785015*^9}, {3.805480860379613*^9, 3.805480975306303*^9}, {
   3.805481055037532*^9, 3.805481138991852*^9}, {3.805481207897601*^9, 
   3.8054813233177137`*^9}, {3.805481668996467*^9, 3.8054816937761745`*^9}, {
   3.8054824596940994`*^9, 3.805482511675572*^9}, 3.80548270058698*^9, {
   3.805482938776911*^9, 3.805482957449237*^9}, 3.8054830084992285`*^9, 
   3.8055023938153105`*^9, 3.805502465059903*^9, 3.8055025620124936`*^9, {
   3.8055025944080496`*^9, 3.805502624323163*^9}, {3.8055026772153587`*^9, 
   3.805502783940276*^9}, {3.805502854239046*^9, 3.80550289141193*^9}, 
   3.805539092142848*^9, 3.8055392008147535`*^9, {3.805539266359166*^9, 
   3.8055392950748816`*^9}, 3.8055476424770594`*^9, {3.805547743548753*^9, 
   3.805547782692649*^9}, 3.806061311121106*^9, {3.8060613762090282`*^9, 
   3.806061399563547*^9}, 3.806061485787599*^9, {3.8060618969884276`*^9, 
   3.8060619413353577`*^9}, 3.8060619883466463`*^9, 3.806062076635193*^9, {
   3.806062121638563*^9, 3.8060622477261877`*^9}, 3.8061436402582865`*^9, 
   3.812479908939375*^9, {3.8124805431737375`*^9, 3.8124805521521053`*^9}, 
   3.8127080260794463`*^9, 3.8127082080786448`*^9, 3.8127082767170734`*^9, 
   3.8127085542792435`*^9, 3.8127103778091154`*^9, 3.81271041590324*^9, 
   3.8127105663518934`*^9, 3.8127106029664683`*^9, 3.8127110945911274`*^9, {
   3.8127116204864073`*^9, 3.8127116438698387`*^9}, 3.8127120512537746`*^9, 
   3.8127132694632845`*^9, 3.8127135282704067`*^9},
 CellLabel->"Out[81]=",ExpressionUUID->"e81b2a24-b6cd-4a79-9cdb-b5f34cfa1aa8"],

Cell[BoxData["\<\"\[CapitalDelta]W\"\>"], "Output",
 CellChangeTimes->{
  3.8054797335618596`*^9, 3.805479768529723*^9, 3.805479962547964*^9, 
   3.805480011263193*^9, 3.805480158628708*^9, 3.8054805207302823`*^9, {
   3.8054806478378954`*^9, 3.8054807227787147`*^9}, {3.805480754776679*^9, 
   3.805480802785015*^9}, {3.805480860379613*^9, 3.805480975306303*^9}, {
   3.805481055037532*^9, 3.805481138991852*^9}, {3.805481207897601*^9, 
   3.8054813233177137`*^9}, {3.805481668996467*^9, 3.8054816937761745`*^9}, {
   3.8054824596940994`*^9, 3.805482511675572*^9}, 3.80548270058698*^9, {
   3.805482938776911*^9, 3.805482957449237*^9}, 3.8054830084992285`*^9, 
   3.8055023938153105`*^9, 3.805502465059903*^9, 3.8055025620124936`*^9, {
   3.8055025944080496`*^9, 3.805502624323163*^9}, {3.8055026772153587`*^9, 
   3.805502783940276*^9}, {3.805502854239046*^9, 3.80550289141193*^9}, 
   3.805539092142848*^9, 3.8055392008147535`*^9, {3.805539266359166*^9, 
   3.8055392950748816`*^9}, 3.8055476424770594`*^9, {3.805547743548753*^9, 
   3.805547782692649*^9}, 3.806061311121106*^9, {3.8060613762090282`*^9, 
   3.806061399563547*^9}, 3.806061485787599*^9, {3.8060618969884276`*^9, 
   3.8060619413353577`*^9}, 3.8060619883466463`*^9, 3.806062076635193*^9, {
   3.806062121638563*^9, 3.8060622477261877`*^9}, 3.8061436402582865`*^9, 
   3.812479908939375*^9, {3.8124805431737375`*^9, 3.8124805521521053`*^9}, 
   3.8127080260794463`*^9, 3.8127082080786448`*^9, 3.8127082767170734`*^9, 
   3.8127085542792435`*^9, 3.8127103778091154`*^9, 3.81271041590324*^9, 
   3.8127105663518934`*^9, 3.8127106029664683`*^9, 3.8127110945911274`*^9, {
   3.8127116204864073`*^9, 3.8127116438698387`*^9}, 3.8127120512537746`*^9, 
   3.8127132694632845`*^9, 3.812713528273399*^9},
 CellLabel->"Out[82]=",ExpressionUUID->"f077f4d8-40d7-4174-9aa6-28f34157ee4a"],

Cell[BoxData["39.0627885503236`"], "Output",
 CellChangeTimes->{
  3.8054797335618596`*^9, 3.805479768529723*^9, 3.805479962547964*^9, 
   3.805480011263193*^9, 3.805480158628708*^9, 3.8054805207302823`*^9, {
   3.8054806478378954`*^9, 3.8054807227787147`*^9}, {3.805480754776679*^9, 
   3.805480802785015*^9}, {3.805480860379613*^9, 3.805480975306303*^9}, {
   3.805481055037532*^9, 3.805481138991852*^9}, {3.805481207897601*^9, 
   3.8054813233177137`*^9}, {3.805481668996467*^9, 3.8054816937761745`*^9}, {
   3.8054824596940994`*^9, 3.805482511675572*^9}, 3.80548270058698*^9, {
   3.805482938776911*^9, 3.805482957449237*^9}, 3.8054830084992285`*^9, 
   3.8055023938153105`*^9, 3.805502465059903*^9, 3.8055025620124936`*^9, {
   3.8055025944080496`*^9, 3.805502624323163*^9}, {3.8055026772153587`*^9, 
   3.805502783940276*^9}, {3.805502854239046*^9, 3.80550289141193*^9}, 
   3.805539092142848*^9, 3.8055392008147535`*^9, {3.805539266359166*^9, 
   3.8055392950748816`*^9}, 3.8055476424770594`*^9, {3.805547743548753*^9, 
   3.805547782692649*^9}, 3.806061311121106*^9, {3.8060613762090282`*^9, 
   3.806061399563547*^9}, 3.806061485787599*^9, {3.8060618969884276`*^9, 
   3.8060619413353577`*^9}, 3.8060619883466463`*^9, 3.806062076635193*^9, {
   3.806062121638563*^9, 3.8060622477261877`*^9}, 3.8061436402582865`*^9, 
   3.812479908939375*^9, {3.8124805431737375`*^9, 3.8124805521521053`*^9}, 
   3.8127080260794463`*^9, 3.8127082080786448`*^9, 3.8127082767170734`*^9, 
   3.8127085542792435`*^9, 3.8127103778091154`*^9, 3.81271041590324*^9, 
   3.8127105663518934`*^9, 3.8127106029664683`*^9, 3.8127110945911274`*^9, {
   3.8127116204864073`*^9, 3.8127116438698387`*^9}, 3.8127120512537746`*^9, 
   3.8127132694632845`*^9, 3.8127135282753935`*^9},
 CellLabel->"Out[83]=",ExpressionUUID->"e7d5684c-c7dc-4ea5-b8df-f6f5987406f1"],

Cell[BoxData["\<\"relative slopes at ehat\"\>"], "Output",
 CellChangeTimes->{
  3.8054797335618596`*^9, 3.805479768529723*^9, 3.805479962547964*^9, 
   3.805480011263193*^9, 3.805480158628708*^9, 3.8054805207302823`*^9, {
   3.8054806478378954`*^9, 3.8054807227787147`*^9}, {3.805480754776679*^9, 
   3.805480802785015*^9}, {3.805480860379613*^9, 3.805480975306303*^9}, {
   3.805481055037532*^9, 3.805481138991852*^9}, {3.805481207897601*^9, 
   3.8054813233177137`*^9}, {3.805481668996467*^9, 3.8054816937761745`*^9}, {
   3.8054824596940994`*^9, 3.805482511675572*^9}, 3.80548270058698*^9, {
   3.805482938776911*^9, 3.805482957449237*^9}, 3.8054830084992285`*^9, 
   3.8055023938153105`*^9, 3.805502465059903*^9, 3.8055025620124936`*^9, {
   3.8055025944080496`*^9, 3.805502624323163*^9}, {3.8055026772153587`*^9, 
   3.805502783940276*^9}, {3.805502854239046*^9, 3.80550289141193*^9}, 
   3.805539092142848*^9, 3.8055392008147535`*^9, {3.805539266359166*^9, 
   3.8055392950748816`*^9}, 3.8055476424770594`*^9, {3.805547743548753*^9, 
   3.805547782692649*^9}, 3.806061311121106*^9, {3.8060613762090282`*^9, 
   3.806061399563547*^9}, 3.806061485787599*^9, {3.8060618969884276`*^9, 
   3.8060619413353577`*^9}, 3.8060619883466463`*^9, 3.806062076635193*^9, {
   3.806062121638563*^9, 3.8060622477261877`*^9}, 3.8061436402582865`*^9, 
   3.812479908939375*^9, {3.8124805431737375`*^9, 3.8124805521521053`*^9}, 
   3.8127080260794463`*^9, 3.8127082080786448`*^9, 3.8127082767170734`*^9, 
   3.8127085542792435`*^9, 3.8127103778091154`*^9, 3.81271041590324*^9, 
   3.8127105663518934`*^9, 3.8127106029664683`*^9, 3.8127110945911274`*^9, {
   3.8127116204864073`*^9, 3.8127116438698387`*^9}, 3.8127120512537746`*^9, 
   3.8127132694632845`*^9, 3.812713528277388*^9},
 CellLabel->"Out[84]=",ExpressionUUID->"e9cc61dc-6752-4f12-80c1-11c137ce3539"],

Cell[BoxData["10.`"], "Output",
 CellChangeTimes->{
  3.8054797335618596`*^9, 3.805479768529723*^9, 3.805479962547964*^9, 
   3.805480011263193*^9, 3.805480158628708*^9, 3.8054805207302823`*^9, {
   3.8054806478378954`*^9, 3.8054807227787147`*^9}, {3.805480754776679*^9, 
   3.805480802785015*^9}, {3.805480860379613*^9, 3.805480975306303*^9}, {
   3.805481055037532*^9, 3.805481138991852*^9}, {3.805481207897601*^9, 
   3.8054813233177137`*^9}, {3.805481668996467*^9, 3.8054816937761745`*^9}, {
   3.8054824596940994`*^9, 3.805482511675572*^9}, 3.80548270058698*^9, {
   3.805482938776911*^9, 3.805482957449237*^9}, 3.8054830084992285`*^9, 
   3.8055023938153105`*^9, 3.805502465059903*^9, 3.8055025620124936`*^9, {
   3.8055025944080496`*^9, 3.805502624323163*^9}, {3.8055026772153587`*^9, 
   3.805502783940276*^9}, {3.805502854239046*^9, 3.80550289141193*^9}, 
   3.805539092142848*^9, 3.8055392008147535`*^9, {3.805539266359166*^9, 
   3.8055392950748816`*^9}, 3.8055476424770594`*^9, {3.805547743548753*^9, 
   3.805547782692649*^9}, 3.806061311121106*^9, {3.8060613762090282`*^9, 
   3.806061399563547*^9}, 3.806061485787599*^9, {3.8060618969884276`*^9, 
   3.8060619413353577`*^9}, 3.8060619883466463`*^9, 3.806062076635193*^9, {
   3.806062121638563*^9, 3.8060622477261877`*^9}, 3.8061436402582865`*^9, 
   3.812479908939375*^9, {3.8124805431737375`*^9, 3.8124805521521053`*^9}, 
   3.8127080260794463`*^9, 3.8127082080786448`*^9, 3.8127082767170734`*^9, 
   3.8127085542792435`*^9, 3.8127103778091154`*^9, 3.81271041590324*^9, 
   3.8127105663518934`*^9, 3.8127106029664683`*^9, 3.8127110945911274`*^9, {
   3.8127116204864073`*^9, 3.8127116438698387`*^9}, 3.8127120512537746`*^9, 
   3.8127132694632845`*^9, 3.8127135282803802`*^9},
 CellLabel->"Out[85]=",ExpressionUUID->"c76bc211-51ba-47d5-8ef6-742f837db588"],

Cell[BoxData["\<\"ESNB True from using Weitzman's tax\"\>"], "Output",
 CellChangeTimes->{
  3.8054797335618596`*^9, 3.805479768529723*^9, 3.805479962547964*^9, 
   3.805480011263193*^9, 3.805480158628708*^9, 3.8054805207302823`*^9, {
   3.8054806478378954`*^9, 3.8054807227787147`*^9}, {3.805480754776679*^9, 
   3.805480802785015*^9}, {3.805480860379613*^9, 3.805480975306303*^9}, {
   3.805481055037532*^9, 3.805481138991852*^9}, {3.805481207897601*^9, 
   3.8054813233177137`*^9}, {3.805481668996467*^9, 3.8054816937761745`*^9}, {
   3.8054824596940994`*^9, 3.805482511675572*^9}, 3.80548270058698*^9, {
   3.805482938776911*^9, 3.805482957449237*^9}, 3.8054830084992285`*^9, 
   3.8055023938153105`*^9, 3.805502465059903*^9, 3.8055025620124936`*^9, {
   3.8055025944080496`*^9, 3.805502624323163*^9}, {3.8055026772153587`*^9, 
   3.805502783940276*^9}, {3.805502854239046*^9, 3.80550289141193*^9}, 
   3.805539092142848*^9, 3.8055392008147535`*^9, {3.805539266359166*^9, 
   3.8055392950748816`*^9}, 3.8055476424770594`*^9, {3.805547743548753*^9, 
   3.805547782692649*^9}, 3.806061311121106*^9, {3.8060613762090282`*^9, 
   3.806061399563547*^9}, 3.806061485787599*^9, {3.8060618969884276`*^9, 
   3.8060619413353577`*^9}, 3.8060619883466463`*^9, 3.806062076635193*^9, {
   3.806062121638563*^9, 3.8060622477261877`*^9}, 3.8061436402582865`*^9, 
   3.812479908939375*^9, {3.8124805431737375`*^9, 3.8124805521521053`*^9}, 
   3.8127080260794463`*^9, 3.8127082080786448`*^9, 3.8127082767170734`*^9, 
   3.8127085542792435`*^9, 3.8127103778091154`*^9, 3.81271041590324*^9, 
   3.8127105663518934`*^9, 3.8127106029664683`*^9, 3.8127110945911274`*^9, {
   3.8127116204864073`*^9, 3.8127116438698387`*^9}, 3.8127120512537746`*^9, 
   3.8127132694632845`*^9, 3.812713528282375*^9},
 CellLabel->"Out[86]=",ExpressionUUID->"2ec0e5ac-177a-4223-9ff7-b82ddf6060a7"],

Cell[BoxData[
 RowBox[{"-", "268.5866667372111`"}]], "Output",
 CellChangeTimes->{
  3.8054797335618596`*^9, 3.805479768529723*^9, 3.805479962547964*^9, 
   3.805480011263193*^9, 3.805480158628708*^9, 3.8054805207302823`*^9, {
   3.8054806478378954`*^9, 3.8054807227787147`*^9}, {3.805480754776679*^9, 
   3.805480802785015*^9}, {3.805480860379613*^9, 3.805480975306303*^9}, {
   3.805481055037532*^9, 3.805481138991852*^9}, {3.805481207897601*^9, 
   3.8054813233177137`*^9}, {3.805481668996467*^9, 3.8054816937761745`*^9}, {
   3.8054824596940994`*^9, 3.805482511675572*^9}, 3.80548270058698*^9, {
   3.805482938776911*^9, 3.805482957449237*^9}, 3.8054830084992285`*^9, 
   3.8055023938153105`*^9, 3.805502465059903*^9, 3.8055025620124936`*^9, {
   3.8055025944080496`*^9, 3.805502624323163*^9}, {3.8055026772153587`*^9, 
   3.805502783940276*^9}, {3.805502854239046*^9, 3.80550289141193*^9}, 
   3.805539092142848*^9, 3.8055392008147535`*^9, {3.805539266359166*^9, 
   3.8055392950748816`*^9}, 3.8055476424770594`*^9, {3.805547743548753*^9, 
   3.805547782692649*^9}, 3.806061311121106*^9, {3.8060613762090282`*^9, 
   3.806061399563547*^9}, 3.806061485787599*^9, {3.8060618969884276`*^9, 
   3.8060619413353577`*^9}, 3.8060619883466463`*^9, 3.806062076635193*^9, {
   3.806062121638563*^9, 3.8060622477261877`*^9}, 3.8061436402582865`*^9, 
   3.812479908939375*^9, {3.8124805431737375`*^9, 3.8124805521521053`*^9}, 
   3.8127080260794463`*^9, 3.8127082080786448`*^9, 3.8127082767170734`*^9, 
   3.8127085542792435`*^9, 3.8127103778091154`*^9, 3.81271041590324*^9, 
   3.8127105663518934`*^9, 3.8127106029664683`*^9, 3.8127110945911274`*^9, {
   3.8127116204864073`*^9, 3.8127116438698387`*^9}, 3.8127120512537746`*^9, 
   3.8127132694632845`*^9, 3.8127135282843695`*^9},
 CellLabel->"Out[87]=",ExpressionUUID->"8d119c0a-9df6-498b-8edd-b29362c39aa0"],

Cell[BoxData["\<\"True Expected Gain in SNB from regulating under Weitzman's \
tax\"\>"], "Output",
 CellChangeTimes->{
  3.8054797335618596`*^9, 3.805479768529723*^9, 3.805479962547964*^9, 
   3.805480011263193*^9, 3.805480158628708*^9, 3.8054805207302823`*^9, {
   3.8054806478378954`*^9, 3.8054807227787147`*^9}, {3.805480754776679*^9, 
   3.805480802785015*^9}, {3.805480860379613*^9, 3.805480975306303*^9}, {
   3.805481055037532*^9, 3.805481138991852*^9}, {3.805481207897601*^9, 
   3.8054813233177137`*^9}, {3.805481668996467*^9, 3.8054816937761745`*^9}, {
   3.8054824596940994`*^9, 3.805482511675572*^9}, 3.80548270058698*^9, {
   3.805482938776911*^9, 3.805482957449237*^9}, 3.8054830084992285`*^9, 
   3.8055023938153105`*^9, 3.805502465059903*^9, 3.8055025620124936`*^9, {
   3.8055025944080496`*^9, 3.805502624323163*^9}, {3.8055026772153587`*^9, 
   3.805502783940276*^9}, {3.805502854239046*^9, 3.80550289141193*^9}, 
   3.805539092142848*^9, 3.8055392008147535`*^9, {3.805539266359166*^9, 
   3.8055392950748816`*^9}, 3.8055476424770594`*^9, {3.805547743548753*^9, 
   3.805547782692649*^9}, 3.806061311121106*^9, {3.8060613762090282`*^9, 
   3.806061399563547*^9}, 3.806061485787599*^9, {3.8060618969884276`*^9, 
   3.8060619413353577`*^9}, 3.8060619883466463`*^9, 3.806062076635193*^9, {
   3.806062121638563*^9, 3.8060622477261877`*^9}, 3.8061436402582865`*^9, 
   3.812479908939375*^9, {3.8124805431737375`*^9, 3.8124805521521053`*^9}, 
   3.8127080260794463`*^9, 3.8127082080786448`*^9, 3.8127082767170734`*^9, 
   3.8127085542792435`*^9, 3.8127103778091154`*^9, 3.81271041590324*^9, 
   3.8127105663518934`*^9, 3.8127106029664683`*^9, 3.8127110945911274`*^9, {
   3.8127116204864073`*^9, 3.8127116438698387`*^9}, 3.8127120512537746`*^9, 
   3.8127132694632845`*^9, 3.8127135282883854`*^9},
 CellLabel->"Out[88]=",ExpressionUUID->"a810af57-ffee-49d3-adf1-99b3dac9e74d"],

Cell[BoxData["15288.310739767232`"], "Output",
 CellChangeTimes->{
  3.8054797335618596`*^9, 3.805479768529723*^9, 3.805479962547964*^9, 
   3.805480011263193*^9, 3.805480158628708*^9, 3.8054805207302823`*^9, {
   3.8054806478378954`*^9, 3.8054807227787147`*^9}, {3.805480754776679*^9, 
   3.805480802785015*^9}, {3.805480860379613*^9, 3.805480975306303*^9}, {
   3.805481055037532*^9, 3.805481138991852*^9}, {3.805481207897601*^9, 
   3.8054813233177137`*^9}, {3.805481668996467*^9, 3.8054816937761745`*^9}, {
   3.8054824596940994`*^9, 3.805482511675572*^9}, 3.80548270058698*^9, {
   3.805482938776911*^9, 3.805482957449237*^9}, 3.8054830084992285`*^9, 
   3.8055023938153105`*^9, 3.805502465059903*^9, 3.8055025620124936`*^9, {
   3.8055025944080496`*^9, 3.805502624323163*^9}, {3.8055026772153587`*^9, 
   3.805502783940276*^9}, {3.805502854239046*^9, 3.80550289141193*^9}, 
   3.805539092142848*^9, 3.8055392008147535`*^9, {3.805539266359166*^9, 
   3.8055392950748816`*^9}, 3.8055476424770594`*^9, {3.805547743548753*^9, 
   3.805547782692649*^9}, 3.806061311121106*^9, {3.8060613762090282`*^9, 
   3.806061399563547*^9}, 3.806061485787599*^9, {3.8060618969884276`*^9, 
   3.8060619413353577`*^9}, 3.8060619883466463`*^9, 3.806062076635193*^9, {
   3.806062121638563*^9, 3.8060622477261877`*^9}, 3.8061436402582865`*^9, 
   3.812479908939375*^9, {3.8124805431737375`*^9, 3.8124805521521053`*^9}, 
   3.8127080260794463`*^9, 3.8127082080786448`*^9, 3.8127082767170734`*^9, 
   3.8127085542792435`*^9, 3.8127103778091154`*^9, 3.81271041590324*^9, 
   3.8127105663518934`*^9, 3.8127106029664683`*^9, 3.8127110945911274`*^9, {
   3.8127116204864073`*^9, 3.8127116438698387`*^9}, 3.8127120512537746`*^9, 
   3.8127132694632845`*^9, 3.81271352829138*^9},
 CellLabel->"Out[89]=",ExpressionUUID->"31490e04-1c81-4bf2-86d7-499e140d5407"],

Cell[BoxData["\<\"% diff in ESNB from policy based on correlated \
uncertainty\"\>"], "Output",
 CellChangeTimes->{
  3.8054797335618596`*^9, 3.805479768529723*^9, 3.805479962547964*^9, 
   3.805480011263193*^9, 3.805480158628708*^9, 3.8054805207302823`*^9, {
   3.8054806478378954`*^9, 3.8054807227787147`*^9}, {3.805480754776679*^9, 
   3.805480802785015*^9}, {3.805480860379613*^9, 3.805480975306303*^9}, {
   3.805481055037532*^9, 3.805481138991852*^9}, {3.805481207897601*^9, 
   3.8054813233177137`*^9}, {3.805481668996467*^9, 3.8054816937761745`*^9}, {
   3.8054824596940994`*^9, 3.805482511675572*^9}, 3.80548270058698*^9, {
   3.805482938776911*^9, 3.805482957449237*^9}, 3.8054830084992285`*^9, 
   3.8055023938153105`*^9, 3.805502465059903*^9, 3.8055025620124936`*^9, {
   3.8055025944080496`*^9, 3.805502624323163*^9}, {3.8055026772153587`*^9, 
   3.805502783940276*^9}, {3.805502854239046*^9, 3.80550289141193*^9}, 
   3.805539092142848*^9, 3.8055392008147535`*^9, {3.805539266359166*^9, 
   3.8055392950748816`*^9}, 3.8055476424770594`*^9, {3.805547743548753*^9, 
   3.805547782692649*^9}, 3.806061311121106*^9, {3.8060613762090282`*^9, 
   3.806061399563547*^9}, 3.806061485787599*^9, {3.8060618969884276`*^9, 
   3.8060619413353577`*^9}, 3.8060619883466463`*^9, 3.806062076635193*^9, {
   3.806062121638563*^9, 3.8060622477261877`*^9}, 3.8061436402582865`*^9, 
   3.812479908939375*^9, {3.8124805431737375`*^9, 3.8124805521521053`*^9}, 
   3.8127080260794463`*^9, 3.8127082080786448`*^9, 3.8127082767170734`*^9, 
   3.8127085542792435`*^9, 3.8127103778091154`*^9, 3.81271041590324*^9, 
   3.8127105663518934`*^9, 3.8127106029664683`*^9, 3.8127110945911274`*^9, {
   3.8127116204864073`*^9, 3.8127116438698387`*^9}, 3.8127120512537746`*^9, 
   3.8127132694632845`*^9, 3.812713528293346*^9},
 CellLabel->"Out[90]=",ExpressionUUID->"c2f1511f-08b1-401a-9607-cee00c5cce0d"],

Cell[BoxData["110.20629333601401`"], "Output",
 CellChangeTimes->{
  3.8054797335618596`*^9, 3.805479768529723*^9, 3.805479962547964*^9, 
   3.805480011263193*^9, 3.805480158628708*^9, 3.8054805207302823`*^9, {
   3.8054806478378954`*^9, 3.8054807227787147`*^9}, {3.805480754776679*^9, 
   3.805480802785015*^9}, {3.805480860379613*^9, 3.805480975306303*^9}, {
   3.805481055037532*^9, 3.805481138991852*^9}, {3.805481207897601*^9, 
   3.8054813233177137`*^9}, {3.805481668996467*^9, 3.8054816937761745`*^9}, {
   3.8054824596940994`*^9, 3.805482511675572*^9}, 3.80548270058698*^9, {
   3.805482938776911*^9, 3.805482957449237*^9}, 3.8054830084992285`*^9, 
   3.8055023938153105`*^9, 3.805502465059903*^9, 3.8055025620124936`*^9, {
   3.8055025944080496`*^9, 3.805502624323163*^9}, {3.8055026772153587`*^9, 
   3.805502783940276*^9}, {3.805502854239046*^9, 3.80550289141193*^9}, 
   3.805539092142848*^9, 3.8055392008147535`*^9, {3.805539266359166*^9, 
   3.8055392950748816`*^9}, 3.8055476424770594`*^9, {3.805547743548753*^9, 
   3.805547782692649*^9}, 3.806061311121106*^9, {3.8060613762090282`*^9, 
   3.806061399563547*^9}, 3.806061485787599*^9, {3.8060618969884276`*^9, 
   3.8060619413353577`*^9}, 3.8060619883466463`*^9, 3.806062076635193*^9, {
   3.806062121638563*^9, 3.8060622477261877`*^9}, 3.8061436402582865`*^9, 
   3.812479908939375*^9, {3.8124805431737375`*^9, 3.8124805521521053`*^9}, 
   3.8127080260794463`*^9, 3.8127082080786448`*^9, 3.8127082767170734`*^9, 
   3.8127085542792435`*^9, 3.8127103778091154`*^9, 3.81271041590324*^9, 
   3.8127105663518934`*^9, 3.8127106029664683`*^9, 3.8127110945911274`*^9, {
   3.8127116204864073`*^9, 3.8127116438698387`*^9}, 3.8127120512537746`*^9, 
   3.8127132694632845`*^9, 3.81271352829534*^9},
 CellLabel->"Out[91]=",ExpressionUUID->"d11a2b31-c7ce-4c9b-ae23-0b248898204c"],

Cell[BoxData["\<\"% diff in gain in ESNB from policy based on correlated \
uncertainty\"\>"], "Output",
 CellChangeTimes->{
  3.8054797335618596`*^9, 3.805479768529723*^9, 3.805479962547964*^9, 
   3.805480011263193*^9, 3.805480158628708*^9, 3.8054805207302823`*^9, {
   3.8054806478378954`*^9, 3.8054807227787147`*^9}, {3.805480754776679*^9, 
   3.805480802785015*^9}, {3.805480860379613*^9, 3.805480975306303*^9}, {
   3.805481055037532*^9, 3.805481138991852*^9}, {3.805481207897601*^9, 
   3.8054813233177137`*^9}, {3.805481668996467*^9, 3.8054816937761745`*^9}, {
   3.8054824596940994`*^9, 3.805482511675572*^9}, 3.80548270058698*^9, {
   3.805482938776911*^9, 3.805482957449237*^9}, 3.8054830084992285`*^9, 
   3.8055023938153105`*^9, 3.805502465059903*^9, 3.8055025620124936`*^9, {
   3.8055025944080496`*^9, 3.805502624323163*^9}, {3.8055026772153587`*^9, 
   3.805502783940276*^9}, {3.805502854239046*^9, 3.80550289141193*^9}, 
   3.805539092142848*^9, 3.8055392008147535`*^9, {3.805539266359166*^9, 
   3.8055392950748816`*^9}, 3.8055476424770594`*^9, {3.805547743548753*^9, 
   3.805547782692649*^9}, 3.806061311121106*^9, {3.8060613762090282`*^9, 
   3.806061399563547*^9}, 3.806061485787599*^9, {3.8060618969884276`*^9, 
   3.8060619413353577`*^9}, 3.8060619883466463`*^9, 3.806062076635193*^9, {
   3.806062121638563*^9, 3.8060622477261877`*^9}, 3.8061436402582865`*^9, 
   3.812479908939375*^9, {3.8124805431737375`*^9, 3.8124805521521053`*^9}, 
   3.8127080260794463`*^9, 3.8127082080786448`*^9, 3.8127082767170734`*^9, 
   3.8127085542792435`*^9, 3.8127103778091154`*^9, 3.81271041590324*^9, 
   3.8127105663518934`*^9, 3.8127106029664683`*^9, 3.8127110945911274`*^9, {
   3.8127116204864073`*^9, 3.8127116438698387`*^9}, 3.8127120512537746`*^9, 
   3.8127132694632845`*^9, 3.8127135282973347`*^9},
 CellLabel->"Out[92]=",ExpressionUUID->"dfdd6a26-1654-4507-b357-415e362c9563"],

Cell[BoxData["15.945071800849888`"], "Output",
 CellChangeTimes->{
  3.8054797335618596`*^9, 3.805479768529723*^9, 3.805479962547964*^9, 
   3.805480011263193*^9, 3.805480158628708*^9, 3.8054805207302823`*^9, {
   3.8054806478378954`*^9, 3.8054807227787147`*^9}, {3.805480754776679*^9, 
   3.805480802785015*^9}, {3.805480860379613*^9, 3.805480975306303*^9}, {
   3.805481055037532*^9, 3.805481138991852*^9}, {3.805481207897601*^9, 
   3.8054813233177137`*^9}, {3.805481668996467*^9, 3.8054816937761745`*^9}, {
   3.8054824596940994`*^9, 3.805482511675572*^9}, 3.80548270058698*^9, {
   3.805482938776911*^9, 3.805482957449237*^9}, 3.8054830084992285`*^9, 
   3.8055023938153105`*^9, 3.805502465059903*^9, 3.8055025620124936`*^9, {
   3.8055025944080496`*^9, 3.805502624323163*^9}, {3.8055026772153587`*^9, 
   3.805502783940276*^9}, {3.805502854239046*^9, 3.80550289141193*^9}, 
   3.805539092142848*^9, 3.8055392008147535`*^9, {3.805539266359166*^9, 
   3.8055392950748816`*^9}, 3.8055476424770594`*^9, {3.805547743548753*^9, 
   3.805547782692649*^9}, 3.806061311121106*^9, {3.8060613762090282`*^9, 
   3.806061399563547*^9}, 3.806061485787599*^9, {3.8060618969884276`*^9, 
   3.8060619413353577`*^9}, 3.8060619883466463`*^9, 3.806062076635193*^9, {
   3.806062121638563*^9, 3.8060622477261877`*^9}, 3.8061436402582865`*^9, 
   3.812479908939375*^9, {3.8124805431737375`*^9, 3.8124805521521053`*^9}, 
   3.8127080260794463`*^9, 3.8127082080786448`*^9, 3.8127082767170734`*^9, 
   3.8127085542792435`*^9, 3.8127103778091154`*^9, 3.81271041590324*^9, 
   3.8127105663518934`*^9, 3.8127106029664683`*^9, 3.8127110945911274`*^9, {
   3.8127116204864073`*^9, 3.8127116438698387`*^9}, 3.8127120512537746`*^9, 
   3.8127132694632845`*^9, 3.8127135283003263`*^9},
 CellLabel->"Out[93]=",ExpressionUUID->"9be7a0bb-5585-49cf-8b72-491017429d71"]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{
  RowBox[{"(*", 
   RowBox[{"More", " ", "numerical", " ", "results"}], "*)"}], 
  "\[IndentingNewLine]", 
  RowBox[{"\"\<Components of \[CapitalDelta]\>\"", "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"\[Sigma]2\[Eta]1", "=", 
     RowBox[{
      RowBox[{
       RowBox[{"(", 
        RowBox[{"\[Eta]1M", "-", "\[Eta]1m"}], ")"}], "^", "2"}], "/", 
      "12"}]}], ";"}], "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"\[Sigma]2\[Alpha]", "=", 
     RowBox[{"(", 
      RowBox[{"\[Sigma]2\[Eta]1", "+", 
       RowBox[{"\[Sigma]2\[Eta]2", "*", 
        RowBox[{
         RowBox[{"(", 
          RowBox[{"\[Pi]e2", "*", "ehat"}], ")"}], "^", "2"}]}]}], ")"}]}], 
    ";"}], "\[IndentingNewLine]", "\[IndentingNewLine]", 
   "\"\<Adjustment to \[CapitalDelta]W\>\"", "\[IndentingNewLine]", 
   RowBox[{"W\[CapitalDelta]adjust", "=", 
    RowBox[{"\[Sigma]2\[Eta]2", "*", "\[Sigma]2\[Alpha]", "*", 
     RowBox[{
      RowBox[{"(", 
       RowBox[{"ESNB\[Tau]W", "-", "ESNBhat"}], ")"}], "/", 
      RowBox[{"(", 
       RowBox[{"2", "*", 
        RowBox[{"\[Pi]e2", "^", "2"}]}], ")"}]}]}]}], "\[IndentingNewLine]", 
   "\[IndentingNewLine]", "\"\<Risk Adjustment\>\"", "\[IndentingNewLine]", 
   RowBox[{"Radjust", "=", 
    RowBox[{
     RowBox[{"-", 
      RowBox[{"(", 
       RowBox[{"2", "*", 
        RowBox[{"d1", "/", 
         RowBox[{"\[Pi]e2", "^", "2"}]}]}], ")"}]}], "*", 
     RowBox[{"(", 
      RowBox[{"\[Sigma]2\[Eta]1", "+", 
       RowBox[{"\[Sigma]2\[Eta]2", "*", 
        RowBox[{
         RowBox[{"(", 
          RowBox[{"\[Pi]e2", "*", "ehat"}], ")"}], "^", "2"}]}]}], ")"}]}]}], 
   "\[IndentingNewLine]", "\[IndentingNewLine]", "\"\<Policy Adjustment\>\"", 
   "\[IndentingNewLine]", 
   RowBox[{"Padjust", "=", 
    RowBox[{
     RowBox[{"-", 
      RowBox[{"(", 
       RowBox[{"1", "+", 
        RowBox[{
         RowBox[{"Log", "[", 
          RowBox[{"\[Eta]2M", "/", "\[Eta]2m"}], "]"}], "*", 
         RowBox[{"(", 
          RowBox[{"1", "/", 
           RowBox[{"(", 
            RowBox[{"\[Eta]2M", "-", "\[Eta]2m"}], ")"}]}], ")"}]}]}], 
       ")"}]}], "*", "d1", "*", "ehat", "*", 
     RowBox[{
      RowBox[{"(", 
       RowBox[{"\[Tau]star", "-", "\[Tau]W"}], ")"}], "/", 
      RowBox[{"(", 
       RowBox[{"2", "*", "\[Pi]e2"}], ")"}]}]}]}], "\[IndentingNewLine]", 
   RowBox[{"(*", 
    RowBox[{"Padjust", "=", 
     RowBox[{
      RowBox[{"(", 
       RowBox[{
        RowBox[{
         RowBox[{"(", 
          RowBox[{"\[Sigma]2\[Eta]1", "+", 
           RowBox[{"\[Pi]e2", "*", "ehat"}]}], ")"}], "*", 
         RowBox[{"(", 
          RowBox[{
           RowBox[{"1", "/", "\[Eta]2m"}], "-", 
           RowBox[{"1", "/", "\[Eta]2M"}]}], ")"}], "*", 
         RowBox[{"(", 
          RowBox[{"1", "/", 
           RowBox[{"(", 
            RowBox[{"\[Eta]2M", "-", "\[Eta]2m"}], ")"}]}], ")"}]}], "+", 
        RowBox[{"\[Pi]e2", "*", "ehat"}], "-", 
        RowBox[{"2", "*", 
         RowBox[{"Log", "[", 
          RowBox[{"\[Eta]2M", "/", "\[Eta]2m"}], "]"}], "*", 
         RowBox[{"(", 
          RowBox[{"1", "/", 
           RowBox[{"(", 
            RowBox[{"\[Eta]2M", "-", "\[Eta]2m"}], ")"}]}], ")"}], "*", 
         "\[Pi]e2", "*", "ehat"}]}], ")"}], "*", "d1", "*", 
      RowBox[{
       RowBox[{"(", 
        RowBox[{"\[Tau]star", "-", "\[Tau]W"}], ")"}], "/", 
       RowBox[{"(", 
        RowBox[{"\[Pi]e2", "^", "2"}], ")"}]}]}]}], "*)"}], 
   "\[IndentingNewLine]", "\[IndentingNewLine]", 
   "\"\<Covariance adjustment\>\"", "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"(", 
     RowBox[{"ESNB\[Tau]", "-", "ESNBhat", "-", 
      RowBox[{"(", 
       RowBox[{"ESNB\[Tau]W", "-", "ESNBhat"}], ")"}]}], ")"}], "-", 
    "W\[CapitalDelta]adjust", "-", "Radjust", "-", "Padjust"}]}]}]], "Input",
 CellChangeTimes->{{3.80614364545619*^9, 3.80614367194915*^9}, {
   3.8061437092559566`*^9, 3.806143709506259*^9}, {3.8061440055617886`*^9, 
   3.806144005621766*^9}, {3.8061440362887745`*^9, 3.8061440689899335`*^9}, {
   3.8061441750414095`*^9, 3.80614427972929*^9}, {3.8061443257255077`*^9, 
   3.8061443503502226`*^9}, {3.806144381058831*^9, 3.806144458507057*^9}, {
   3.806144499063833*^9, 3.8061445388783336`*^9}, {3.8061459583359222`*^9, 
   3.8061459667610106`*^9}, 3.8061459972033873`*^9, {3.806147006884633*^9, 
   3.8061470139116764`*^9}, {3.8061471369998045`*^9, 3.8061471603420086`*^9}, 
   3.806147218005783*^9, {3.806148348262578*^9, 3.8061485892288113`*^9}, {
   3.806148699183951*^9, 3.806148704521948*^9}, {3.806148760801255*^9, 
   3.8061487698789587`*^9}, {3.806148903315751*^9, 3.8061489065975275`*^9}, {
   3.8061489815751877`*^9, 3.806148981764081*^9}, {3.8061507571492395`*^9, 
   3.8061507657328095`*^9}, {3.806151034193465*^9, 3.8061511002248435`*^9}, 
   3.806151278968881*^9, {3.8061514532973876`*^9, 3.806151454401266*^9}, {
   3.812479940800783*^9, 3.8124799474137163`*^9}},
 CellLabel->"In[94]:=",ExpressionUUID->"c681558f-2f6c-46a3-929e-fc63fa8b656e"],

Cell[BoxData["\<\"Components of \[CapitalDelta]\"\>"], "Output",
 CellChangeTimes->{
  3.8061444617557545`*^9, {3.806144504024417*^9, 3.8061445160827*^9}, 
   3.8061460019680233`*^9, 3.8061470146301365`*^9, 3.806147165166222*^9, 
   3.8061472185601325`*^9, {3.8061485771377153`*^9, 3.8061485913168364`*^9}, 
   3.8061487106915474`*^9, 3.8061487780953236`*^9, 3.8061489077667394`*^9, 
   3.806148982534346*^9, 3.8061511009732385`*^9, 3.8061512797520437`*^9, 
   3.8061514551326528`*^9, 3.8124799245051765`*^9, 3.812708028977728*^9, 
   3.812716495381525*^9},
 CellLabel->"Out[94]=",ExpressionUUID->"cffe407a-3fb2-4674-9cf2-1d3e7bbb47fe"],

Cell[BoxData["\<\"Adjustment to \[CapitalDelta]W\"\>"], "Output",
 CellChangeTimes->{
  3.8061444617557545`*^9, {3.806144504024417*^9, 3.8061445160827*^9}, 
   3.8061460019680233`*^9, 3.8061470146301365`*^9, 3.806147165166222*^9, 
   3.8061472185601325`*^9, {3.8061485771377153`*^9, 3.8061485913168364`*^9}, 
   3.8061487106915474`*^9, 3.8061487780953236`*^9, 3.8061489077667394`*^9, 
   3.806148982534346*^9, 3.8061511009732385`*^9, 3.8061512797520437`*^9, 
   3.8061514551326528`*^9, 3.8124799245051765`*^9, 3.812708028977728*^9, 
   3.8127164953895288`*^9},
 CellLabel->"Out[97]=",ExpressionUUID->"9317045c-939d-432d-9b1a-4c8b3dd9c054"],

Cell[BoxData["41.19933913183643`"], "Output",
 CellChangeTimes->{
  3.8061444617557545`*^9, {3.806144504024417*^9, 3.8061445160827*^9}, 
   3.8061460019680233`*^9, 3.8061470146301365`*^9, 3.806147165166222*^9, 
   3.8061472185601325`*^9, {3.8061485771377153`*^9, 3.8061485913168364`*^9}, 
   3.8061487106915474`*^9, 3.8061487780953236`*^9, 3.8061489077667394`*^9, 
   3.806148982534346*^9, 3.8061511009732385`*^9, 3.8061512797520437`*^9, 
   3.8061514551326528`*^9, 3.8124799245051765`*^9, 3.812708028977728*^9, 
   3.812716495392557*^9},
 CellLabel->"Out[98]=",ExpressionUUID->"f1026d01-d01d-4318-80e3-6990e29bc3f9"],

Cell[BoxData["\<\"Risk Adjustment\"\>"], "Output",
 CellChangeTimes->{
  3.8061444617557545`*^9, {3.806144504024417*^9, 3.8061445160827*^9}, 
   3.8061460019680233`*^9, 3.8061470146301365`*^9, 3.806147165166222*^9, 
   3.8061472185601325`*^9, {3.8061485771377153`*^9, 3.8061485913168364`*^9}, 
   3.8061487106915474`*^9, 3.8061487780953236`*^9, 3.8061489077667394`*^9, 
   3.806148982534346*^9, 3.8061511009732385`*^9, 3.8061512797520437`*^9, 
   3.8061514551326528`*^9, 3.8124799245051765`*^9, 3.812708028977728*^9, 
   3.8127164953945484`*^9},
 CellLabel->"Out[99]=",ExpressionUUID->"e557830e-eb9b-4626-9eed-8d3479a075b3"],

Cell[BoxData[
 RowBox[{"-", "703.1301939058169`"}]], "Output",
 CellChangeTimes->{
  3.8061444617557545`*^9, {3.806144504024417*^9, 3.8061445160827*^9}, 
   3.8061460019680233`*^9, 3.8061470146301365`*^9, 3.806147165166222*^9, 
   3.8061472185601325`*^9, {3.8061485771377153`*^9, 3.8061485913168364`*^9}, 
   3.8061487106915474`*^9, 3.8061487780953236`*^9, 3.8061489077667394`*^9, 
   3.806148982534346*^9, 3.8061511009732385`*^9, 3.8061512797520437`*^9, 
   3.8061514551326528`*^9, 3.8124799245051765`*^9, 3.812708028977728*^9, 
   3.8127164953965454`*^9},
 CellLabel->
  "Out[100]=",ExpressionUUID->"f9d52dc9-b208-475e-b2c0-6f9b369b48a7"],

Cell[BoxData["\<\"Policy Adjustment\"\>"], "Output",
 CellChangeTimes->{
  3.8061444617557545`*^9, {3.806144504024417*^9, 3.8061445160827*^9}, 
   3.8061460019680233`*^9, 3.8061470146301365`*^9, 3.806147165166222*^9, 
   3.8061472185601325`*^9, {3.8061485771377153`*^9, 3.8061485913168364`*^9}, 
   3.8061487106915474`*^9, 3.8061487780953236`*^9, 3.8061489077667394`*^9, 
   3.806148982534346*^9, 3.8061511009732385`*^9, 3.8061512797520437`*^9, 
   3.8061514551326528`*^9, 3.8124799245051765`*^9, 3.812708028977728*^9, 
   3.8127164953985095`*^9},
 CellLabel->
  "Out[101]=",ExpressionUUID->"ec7766e0-74e4-4f44-b33d-d57c933d8c9a"],

Cell[BoxData["841.5911720038362`"], "Output",
 CellChangeTimes->{
  3.8061444617557545`*^9, {3.806144504024417*^9, 3.8061445160827*^9}, 
   3.8061460019680233`*^9, 3.8061470146301365`*^9, 3.806147165166222*^9, 
   3.8061472185601325`*^9, {3.8061485771377153`*^9, 3.8061485913168364`*^9}, 
   3.8061487106915474`*^9, 3.8061487780953236`*^9, 3.8061489077667394`*^9, 
   3.806148982534346*^9, 3.8061511009732385`*^9, 3.8061512797520437`*^9, 
   3.8061514551326528`*^9, 3.8124799245051765`*^9, 3.812708028977728*^9, 
   3.8127164954021206`*^9},
 CellLabel->
  "Out[102]=",ExpressionUUID->"2e30abc2-d077-47e7-b83c-728978a40f20"],

Cell[BoxData["\<\"Covariance adjustment\"\>"], "Output",
 CellChangeTimes->{
  3.8061444617557545`*^9, {3.806144504024417*^9, 3.8061445160827*^9}, 
   3.8061460019680233`*^9, 3.8061470146301365`*^9, 3.806147165166222*^9, 
   3.8061472185601325`*^9, {3.8061485771377153`*^9, 3.8061485913168364`*^9}, 
   3.8061487106915474`*^9, 3.8061487780953236`*^9, 3.8061489077667394`*^9, 
   3.806148982534346*^9, 3.8061511009732385`*^9, 3.8061512797520437`*^9, 
   3.8061514551326528`*^9, 3.8124799245051765`*^9, 3.812708028977728*^9, 
   3.812716495405128*^9},
 CellLabel->
  "Out[103]=",ExpressionUUID->"5ae0f7f7-a0a0-4e5f-836f-db7cdb40c732"],

Cell[BoxData[
 RowBox[{"-", "802.5115167265739`"}]], "Output",
 CellChangeTimes->{
  3.8061444617557545`*^9, {3.806144504024417*^9, 3.8061445160827*^9}, 
   3.8061460019680233`*^9, 3.8061470146301365`*^9, 3.806147165166222*^9, 
   3.8061472185601325`*^9, {3.8061485771377153`*^9, 3.8061485913168364`*^9}, 
   3.8061487106915474`*^9, 3.8061487780953236`*^9, 3.8061489077667394`*^9, 
   3.806148982534346*^9, 3.8061511009732385`*^9, 3.8061512797520437`*^9, 
   3.8061514551326528`*^9, 3.8124799245051765`*^9, 3.812708028977728*^9, 
   3.812716495407133*^9},
 CellLabel->
  "Out[104]=",ExpressionUUID->"a89aa0ce-7af7-4254-a393-e92d3e676c92"]
}, Open  ]]
},
WindowSize->{1280, 941},
WindowMargins->{{-8, Automatic}, {Automatic, -8}},
FrontEndVersion->"12.0 for Microsoft Windows (64-bit) (April 8, 2019)",
StyleDefinitions->"Default.nb"
]
(* End of Notebook Content *)

(* Internal cache information *)
(*CellTagsOutline
CellTagsIndex->{}
*)
(*CellTagsIndex
CellTagsIndex->{}
*)
(*NotebookFileOutline
Notebook[{
Cell[CellGroupData[{
Cell[580, 22, 945, 21, 67, "Input",ExpressionUUID->"53974844-04cd-41a6-a1d0-a926fb069df5"],
Cell[1528, 45, 843, 12, 32, "Output",ExpressionUUID->"1e6e590f-bfc1-41dc-9797-22260b874c3d"]
}, Open  ]],
Cell[2386, 60, 1746, 37, 86, "Input",ExpressionUUID->"1fc03bfa-a798-4792-97da-c543b3a20ed5"],
Cell[4135, 99, 4049, 95, 257, "Input",ExpressionUUID->"cdf1330f-e9ba-45f1-8a3e-1254eb87e3c4"],
Cell[CellGroupData[{
Cell[8209, 198, 16509, 412, 1131, "Input",ExpressionUUID->"3a775e5e-012d-4f93-9d76-951eebbda49e"],
Cell[24721, 612, 539, 11, 21, "Message",ExpressionUUID->"ff7be5fb-cd86-40ff-829b-8fa8e66dc96e"]
}, Open  ]],
Cell[25275, 626, 6118, 154, 390, "Input",ExpressionUUID->"701aced9-601e-49c4-9f34-1b9fcca91619"],
Cell[31396, 782, 1577, 37, 67, "Input",ExpressionUUID->"eded3f86-e60f-4f25-a0af-902672acae05"],
Cell[32976, 821, 967, 26, 181, "Input",ExpressionUUID->"e26b58c6-972a-48b9-8b6d-3ed929996a4b"],
Cell[CellGroupData[{
Cell[33968, 851, 7877, 154, 1283, "Input",ExpressionUUID->"daa7a5ae-3fc1-4c4c-b9ab-118cd11e24be"],
Cell[41848, 1007, 1842, 25, 32, "Output",ExpressionUUID->"a027acb6-16ca-40c1-8a52-5f3b40c773be"],
Cell[43693, 1034, 1851, 26, 32, "Output",ExpressionUUID->"9d51ae94-dd48-4adb-b153-8c131a503335"],
Cell[45547, 1062, 1887, 26, 32, "Output",ExpressionUUID->"fba2204a-19a0-47ab-aba7-be81f1c8ade0"],
Cell[47437, 1090, 1836, 25, 32, "Output",ExpressionUUID->"f7094cb9-7693-4d26-9455-3d78a1757817"],
Cell[49276, 1117, 1834, 25, 32, "Output",ExpressionUUID->"62f62a72-77ff-4311-bcc5-6bec6e75c0f2"],
Cell[51113, 1144, 1834, 25, 32, "Output",ExpressionUUID->"4c3d9b50-8f79-4e1e-9069-6b9d580a9863"],
Cell[52950, 1171, 1883, 26, 32, "Output",ExpressionUUID->"340e7b4c-8838-436f-861c-9d191d0a75ae"],
Cell[54836, 1199, 1837, 25, 32, "Output",ExpressionUUID->"c1a82b20-1a94-423e-85bb-9457ee3bb4d2"],
Cell[56676, 1226, 1864, 25, 32, "Output",ExpressionUUID->"bf40d956-5cfa-42ba-bb08-50a28e3c4a83"],
Cell[58543, 1253, 1852, 26, 32, "Output",ExpressionUUID->"9708c0bb-d6ce-4ee2-bac0-0cbf909b0013"],
Cell[60398, 1281, 1834, 25, 32, "Output",ExpressionUUID->"720e5528-627c-46be-bdb6-ff14cec00545"],
Cell[62235, 1308, 1836, 25, 32, "Output",ExpressionUUID->"1472e81b-1498-4d04-bfbb-a33feedcf84a"],
Cell[64074, 1335, 1825, 25, 32, "Output",ExpressionUUID->"c93e01df-33e2-40c3-8193-a3bc52f2831b"],
Cell[65902, 1362, 1837, 25, 32, "Output",ExpressionUUID->"1877f085-bd5b-4a44-869f-43c27c956ffe"],
Cell[67742, 1389, 1842, 25, 32, "Output",ExpressionUUID->"968c2c9c-142f-4115-8c3f-423088b70b44"],
Cell[69587, 1416, 1837, 25, 32, "Output",ExpressionUUID->"b16690f5-ed57-4a36-9dd8-e24347c95e7f"],
Cell[71427, 1443, 1883, 26, 32, "Output",ExpressionUUID->"e6cd8845-c0a8-48b3-a185-cb698677d5b8"],
Cell[73313, 1471, 1836, 25, 32, "Output",ExpressionUUID->"6f901064-f319-4180-939b-28c74dee252f"],
Cell[75152, 1498, 1846, 25, 32, "Output",ExpressionUUID->"db9d170e-6b90-44d3-a61c-a2db2bd71a46"],
Cell[77001, 1525, 1853, 26, 32, "Output",ExpressionUUID->"370ed5ee-b937-4bfa-8195-8d310424a36a"],
Cell[78857, 1553, 1885, 26, 32, "Output",ExpressionUUID->"9b0d1a96-2b64-47a7-9db1-ab7c11bdda8c"],
Cell[80745, 1581, 1836, 25, 32, "Output",ExpressionUUID->"ab5a2342-9bed-4ed7-9b85-c19365b0e2cf"],
Cell[82584, 1608, 1853, 25, 32, "Output",ExpressionUUID->"9e13bdc1-29df-489e-8d09-ae1b3cf0d21b"],
Cell[84440, 1635, 1837, 25, 32, "Output",ExpressionUUID->"122c57ef-ccb7-44e4-896a-f9cd3d9dd7aa"],
Cell[86280, 1662, 1840, 25, 32, "Output",ExpressionUUID->"db7abda0-f71c-41a2-9a62-ab41ce886a5e"],
Cell[88123, 1689, 1835, 25, 32, "Output",ExpressionUUID->"df710ae5-1710-4816-b52b-c1771bb5724f"],
Cell[89961, 1716, 1938, 26, 32, "Output",ExpressionUUID->"73d9a939-b8b8-4af2-be6b-ca643ccbdf2e"],
Cell[91902, 1744, 1836, 25, 32, "Output",ExpressionUUID->"e81b2a24-b6cd-4a79-9cdb-b5f34cfa1aa8"],
Cell[93741, 1771, 1840, 25, 32, "Output",ExpressionUUID->"f077f4d8-40d7-4174-9aa6-28f34157ee4a"],
Cell[95584, 1798, 1835, 25, 32, "Output",ExpressionUUID->"e7d5684c-c7dc-4ea5-b8df-f6f5987406f1"],
Cell[97422, 1825, 1847, 25, 32, "Output",ExpressionUUID->"e9cc61dc-6752-4f12-80c1-11c137ce3539"],
Cell[99272, 1852, 1822, 25, 32, "Output",ExpressionUUID->"c76bc211-51ba-47d5-8ef6-742f837db588"],
Cell[101097, 1879, 1859, 25, 32, "Output",ExpressionUUID->"2ec0e5ac-177a-4223-9ff7-b82ddf6060a7"],
Cell[102959, 1906, 1853, 26, 32, "Output",ExpressionUUID->"8d119c0a-9df6-498b-8edd-b29362c39aa0"],
Cell[104815, 1934, 1890, 26, 32, "Output",ExpressionUUID->"a810af57-ffee-49d3-adf1-99b3dac9e74d"],
Cell[106708, 1962, 1834, 25, 32, "Output",ExpressionUUID->"31490e04-1c81-4bf2-86d7-499e140d5407"],
Cell[108545, 1989, 1884, 26, 32, "Output",ExpressionUUID->"c2f1511f-08b1-401a-9607-cee00c5cce0d"],
Cell[110432, 2017, 1834, 25, 32, "Output",ExpressionUUID->"d11a2b31-c7ce-4c9b-ae23-0b248898204c"],
Cell[112269, 2044, 1894, 26, 32, "Output",ExpressionUUID->"dfdd6a26-1654-4507-b357-415e362c9563"],
Cell[114166, 2072, 1837, 25, 32, "Output",ExpressionUUID->"9be7a0bb-5585-49cf-8b72-491017429d71"]
}, Open  ]],
Cell[CellGroupData[{
Cell[116040, 2102, 4979, 119, 333, "Input",ExpressionUUID->"c681558f-2f6c-46a3-929e-fc63fa8b656e"],
Cell[121022, 2223, 636, 9, 32, "Output",ExpressionUUID->"cffe407a-3fb2-4674-9cf2-1d3e7bbb47fe"],
Cell[121661, 2234, 639, 9, 32, "Output",ExpressionUUID->"9317045c-939d-432d-9b1a-4c8b3dd9c054"],
Cell[122303, 2245, 617, 9, 32, "Output",ExpressionUUID->"f1026d01-d01d-4318-80e3-6990e29bc3f9"],
Cell[122923, 2256, 624, 9, 32, "Output",ExpressionUUID->"e557830e-eb9b-4626-9eed-8d3479a075b3"],
Cell[123550, 2267, 640, 11, 32, "Output",ExpressionUUID->"f9d52dc9-b208-475e-b2c0-6f9b369b48a7"],
Cell[124193, 2280, 630, 10, 32, "Output",ExpressionUUID->"ec7766e0-74e4-4f44-b33d-d57c933d8c9a"],
Cell[124826, 2292, 623, 10, 32, "Output",ExpressionUUID->"2e30abc2-d077-47e7-b83c-728978a40f20"],
Cell[125452, 2304, 632, 10, 32, "Output",ExpressionUUID->"5ae0f7f7-a0a0-4e5f-836f-db7cdb40c732"],
Cell[126087, 2316, 638, 11, 32, "Output",ExpressionUUID->"a89aa0ce-7af7-4254-a393-e92d3e676c92"]
}, Open  ]]
}
]
*)

